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1 Introduction 


pyFEMM is a Python package that allows for the operation of Finite Element Method Magnetics 
(FEMM) via a library of Python functions. The module only works on Windows because it relies on the 
win32com module for communication with FEMM via FEMM’s out-of-proc ActiveX server. 


When pyFEMM starts up a FEMM process, the usual FEMM user interface is displayed and is fully 
functional. The user then has the choice of accomplishing modeling and analysis tasks either 
exclusively through functions implemented by the toolbox, or by a combination of manual and 
programmatic operations — whichever is easiest for the task at hand. 


This document contains a detailed description of all pyFEMM commands. The syntax of the PWFEMM 
toolbox closely mirrors that of FEMM’s existing Lua and Matlab/Octave scripting interfaces. 


2 Installation and Startup 

2.1 Installation / Dependencies 

The package can be installed via pip: 
pip install pyfemm 


from the command line. The package also needs the win32com model (part of pywin32) to 
communicate between Python and FEMM. The pywin32 package should be installed automatically 
as a dependency of pyFEMM. However, if manual installation of the package is desired, it can also 
be done via pip from the command line: 


pip install pypiwin32 


Since the Matplotlib plotting package is used in some pyFEMM examples, it is also listed as a 
dependency and installed with pyFEMM. Again, if manual installation of Matplotlib is desired, it 
can also be don via pip from the command line: 


pip install matplotlib 


2.2 Startup 
To load the pyFEMM module, use: 


>>> import femm 

either at the beginning of a script or from the Python interpreter. 
Then, use the openfemm() function to connect to FEMM: 

>>> femm.openfemm() 


This function starts up a FEMM process to which PyFEMM calls are sent. When done with pyFEMM, 
the FEMM process can be shut down via the femm.closefemm() function. The FEMM process is also 
automatically shut down when Python is closed. 


Note that in all cases, the FEMM-related command must be preceded by the femm. prefix. 


Several examples that use pyFEMM to analyze various problems are available at the pyFEMM website 
www.femm.info/wiki/pyFEMM. 


3 Common Command Set 


There are several FEMM-specific Octave that are not associated with any particular problem type. 
These functions manipulate the appearance of the main window and other top-level components 
like the Lua console and Point Properties output window. 


e openfemm(bHide) The bHide parameter determines whether or not the main FEMM window is 
visible when FEMM is started. In many cases, it is useful to use openfemm(1) to hide the main 
window while a script is running. In hidden window mode, analysis calls like mi_analyze hide the 
analysis window by default. 


e clearconsole() Clears the output window of the Lua console. 


e newdocument(doctype) Creates a new preprocessor document and opens up a new preprocessor 
window. Specify doctype to be 0 for a magnetics problem, 1 for an electrostatics problem, 2 for a 
heat flow problem, or 3 for a current flow problem. Alternative syntax for this function is 
create(doctype) 


e hideconsole() Hides the floating Lua console window. 


e Hidepointprops() Hides the floating FEMM Properties display window. 
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messagebox(’message’) displays the ’message’ string to the screen in a pop-up message box. 
opendocument(’filename’) Opens a document specified by filename. 


print(item1,item2,...) This is standard Lua “print” command directed to the output of the Lua 
console window in FEMM. Any number of comma-separated items can be printed at once via the 
print command. 


prompt(’message’) This function allows a FEMM script to prompt a user for input. When this 
command is used, a dialog box pops up with the ’message’ string on the title bar of the dialog box. 
The user can enter in a single line of input via the dialog box. Prompt returns the user’s input to 
Octave and parses it using Octave’s eval command. 


showconsole() Displays the floating Lua console window. 
showpointprops() Displays the floating FEMM Properties display window. 


smartmesh(state) Calling with a state of O turns off “smart mesh” functionality for the present 
session; calling with a state of 1 turns “smarth meshing” on. The smartmesh function applies across 
all problem types. To set smart meshing on or off on a file-by-file basis, use 
mi_smartmesh,ei_smartmesh, hi_smartmesh, and ci_smartmesh, depending on problem type. 
Note that calling the global smartmesh with a value of -1 sets the program to defer to the file-by- 
file setting rather than forcing all smart meshing on or off. 


main_minimize() minimizes the main FEMM window. 
main_maximize() maximizes the main FEMM window. 
main_restore() restores the main FEMM window from a minimized or maximized state. 


main_resize(width,height) resizes the main FEMM window client area to width x height. 


Magnetics Preprocessor Command Set 


Several different commands are available in the preprocessor. 


4.1 


Object Add/Remove Commands 


mi_addnode(x,y) Add a new node at x,y 


mi_addsegment(x1,y1,x2,y2) Add a new line segment from node closest to (x1,y1) to node closest 
to (x2,y2) 


mi_addblocklabel(x,y) Add a new block label at (x,y) 
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mi_addarc(x1,y1,x2,y2,angle, maxseg) Add a new arc segment from the nearest node to (x1,y1) to 
the nearest node to (x2,y2) with angle ‘angle’ divided into ‘maxseg’ segments. 


mi_drawline(x1,y1,x2,y2) Adds nodes at (x1,y1) and (x2,y2) and adds a line between the nodes. 


mi_drawarc(x1,y1,x2,y2,angle, maxseg) Adds nodes at (x1,y1) and (x2,y2) and adds an arc of the 
specified angle and discretization connecting the nodes. 


mi_drawrectangle(x1,y1,x2,y2) Adds nodes at the corners of a rectangle defined by the points 
(x1,y1) and (x2,y2), then adds segments connecting the corners of the rectangle. 


mi_deleteselected Delete all selected objects. 
mi_deleteselectednodes Delete selected nodes. 
mi_deleteselectedlabels Delete selected block labels. 
mi_deleteselectedsegments Delete selected segments. 


mi_deleteselectedarcsegments Delete selects arcs. 


Geometry Selection Commands 
mi_clearselected Clear all selected nodes, blocks, segments and arc segments. 
mi_selectsegment(x,y) Select the line segment closest to (x,y) 
mi_selectnode(x,y) Select the node closest to (x,y). Returns the coordinates of the selected node. 
mi_selectlabel(x,y) Select the label closet to (x,y). Returns the coordinates of the selected label. 


mi_selectarcsegment(x,y) Select the arc segment closest to (x,y) 


mi_selectgroup(n) Select the n" group of nodes, segments, arc segments and blocklabels. 
This function will clear all previously selected elements and leave the editmode in 4 (group) 


mi_selectcircle(x,y,R,editmode) selects objects within a circle of radius R centered at (x,y). If only 
x, y, and R paramters are given, the current edit mode is used. If the editmode parameter is used, 
0 denotes nodes, 2 denotes block labels, 2 denotes segments, 3 denotes arcs, and 4 specifies that 
all entity types are to be selected. 


mi_selectrectangle(x1,y1,x2,y2,editmode) selects objects within a rectangle defined by points 
(x1,y1) and (x2,y2). If no editmode parameter is supplied, the current edit mode is used. If the 
editmode parameter is used, 0 denotes nodes, 2 denotes block labels, 2 denotes segments, 3 
denotes arcs, and 4 specifies that all entity types are to be selected. 


4.3 Object Labeling Commands 


e mi_setnodeprop(’propname’,groupno) Set the selected nodes to have the nodal property 
"propname’ and group number groupno. 


e mi_setblockprop(’blockname’, automesh, meshsize, ’incircuit’, magdir, group, turns) Set the 
selected block labels to have the properties: 


Block property ’blockname’. 


— automesh: O = mesher defers to mesh size constraint defined in meshsize, 1 = mesher 
automatically chooses the mesh density. 


meshsize: size constraint on the mesh in the block marked by this label. 


Block is a member of the circuit named ’incircuit’ 


The magnetization is directed along an angle in measured in degrees denoted by the 
parameter magdir 


— A member of group number group 


The number of turns associated with this label is denoted by turns. 


e mi_setsegmentprop(’propname’, elementsize, automesh, hide, group) Set the selected segments 
to have: 


Boundary property ’propname’ 


Local element size along segment no greater than elementsize 


automesh: 0 = mesher defers to the element constraint defined by elementsize, 1 = mesher 
automatically chooses mesh size along the selected segments 


— hide: 0 = not hidden in post-processor, 1 == hidden in post processor 


— A member of group number group 


e mi_setarcsegmentprop(maxsegdeg, ’propname’, hide, group) Set the selected arc segments to: 


Meshed with elements that span at most maxsegdeg degrees per element 


Boundary property ’propname’ 
— hide: 0 = not hidden in post-processor, 1 == hidden in post processor 


— A member of group number group 


e mi_setgroup(n) Set the group associated of the selected items to n. 


4.4 


Problem Commands 


e mi_probdef(freq,units,type,precision,depth,minangle,(acsolver)) changes the problem definition. 
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Set freq to the desired frequency in Hertz. The units parameter specifies the units used for 
measuring length in the problem domain. Valid ‘units’ entries are ’inches’, ‘millimeters’, 
‘centimeters’, ’mils’, ’meters’, and ’micrometers’. Set the parameter problemtype to ’planar’ for 
a 2-D planar problem, or to ’axi’ for an axisymmetric problem. The precision parameter dictates 
the precision required by the solver. For example, entering 1E-8 requires the RMS of the residual 
to be less than 10-8. A fifth parameter, representing the depth of the problem in the into-the- 
page direction for 2-D planar problems. Specify the depth to be zero for axisymmetric problems. 
The sixth parameter represents the minimum angle constraint sent to the mesh generator — 30 
degrees is the usual choice for this parameter. The acsolver parameter specifies which solver is 
to be used for AC problems: 0 for successive approximation, 1 for Newton. 


mi_analyze(flag) runs the magnetics solver. The flag parameter controls whether the solver 
window is visible or minimized. For a visible window, specify 0. For a minimized window, flag 
should be set to 1. If no value is specified for flag, the visibility of the solver is inherited from the 
main window, i.e. if the main window is minimized, the solver runs minimized, too. 


mi_loadsolution loads and displays the solution corresponding to the current geometry. 


mi_setfocus(’documentname’) Switches the magnetics input file upon which commands are to 
act. If more than one magnetics input file is being edited at a time, this command can be used to 
switch between files so that the mutiple files can be operated upon programmatically. 
’documentname’ should contain the name of the desired document as it appears on the window’s 
title bar. 


mi_saveas(’filename’) saves the file with name ‘filename’. 


Mesh Commands 


mi_createmesh runs triangle to create a mesh. Note that this is not a necessary precursor of 
performing an analysis, as mi_analyze will make sure the mesh is up to date before running an 
analysis. The number of elements in the mesh is pushed back onto the lua stack. 


mi_showmesh shows the mesh. 


mi_purgemesh clears the mesh out of both the screen and memory. 


Editing Commands 


mi_copyrotate(bx, by, angle, copies ) 


— bx, by — base point for rotation 


— angle — angle by which the selected objects are incrementally shifted to make each copy. 
angle is measured in degrees. 


— copies — number of copies to be produced from the selected objects. 
e mi_copyrotate2(bx, by, angle, copies, editaction ) 


— bx, by — base point for rotation 


angle — angle by which the selected objects are incrementally shifted to make each copy. 
angle is measured in degrees. 


— copies — number of copies to be produced from the selected objects. 


editaction 0 —nodes, 1 — lines (segments), 2 —block labels, 3 — arc segments, 4- group 
e mi_copytranslate(dx, dy, copies) 


— dx,dy — distance by which the selected objects are incrementally shifted. 


— copies — number of copies to be produced from the selected objects. 
e mi_copytranslate2(dx, dy, copies, editaction) 


— dx,dy — distance by which the selected objects are incrementally shifted. 


— copies — number of copies to be produced from the selected objects. 
editaction 0 —nodes, 1 — lines (segments), 2 —block labels, 3 — arc segments, 4- group 
e mi_createradius(x,y,r) turns a corner located at (x,y) into a curve of radius r. 
e mi_moverotate(bx,by,shiftangle) 


— bx, by — base point for rotation 


— shiftangle — angle in degrees by which the selected objects are rotated. 
e mi_moverotate2(bx,by,shiftangle, editaction) 


— bx, by — base point for rotation 
— shiftangle — angle in degrees by which the selected objects are rotated. 


— editaction 0 —-nodes, 1 — lines (segments), 2 —block labels, 3 — arc segments, 4- group 
e mi_movetranslate(dx,dy) 
— dx,dy — distance by which the selected objects are shifted. 


e mi_movetranslate2(dx,dy,editaction) 
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— dx,dy — distance by which the selected objects are shifted. 


— editaction 0 —-nodes, 1 — lines (segments), 2 —block labels, 3 — arc segments, 4- group 
mi_scale(bx, by,scalefactor) 

— bx, by — base point for scaling 

— scalefactor — a multiplier that determines how much the selected objects are scaled 
mi_scale2(bx,by,scalefactor,editaction) 

— bx, by — base point for scaling 


— scalefactor — a multiplier that determines how much the selected objects are scaled 


— editaction 0 —nodes, 1 — lines (segments), 2 —block labels, 3 — arc segments, 4- group 


mi_mirror(x1,y1,x2,y2) mirror the selected objects about a line passing through the points (x1,y1) 
and (x2,y2). 


mi_mirror2(x1,y1,x2,y2,editaction) mirror the selected objects about a line passing through the 
points (x1,y1) and (x2,y2). Valid editaction entries are O for nodes, 1 for lines (segments), 2 for 
block labels, 3 for arc segments, and 4 for groups. 


mi_seteditmode(editmode) Sets the current editmode to: 


— ’nodes’ - nodes 
— ’segments’ - line segments 


— ’arcsegments’ - arc segments 
— ’blocks’ - block labels 


— ’group’ - selected group 


This command will affect all subsequent uses of the other editing commands if they are used 
without the editaction parameter. 


Zoom Commands 


|” 


mi_zoomnatural zooms to a “natural” view with sensible extents. 
mi_zoomout zooms out by a factor of 50%. 


mi_zoomin zoom in by a factor of 200%. 


4.8 
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mi_zoom(x1,y1,x2,y2) Set the display area to be from the bottom left corner specified by (x1,y1) 
to the top right corner specified by (x2,y2). 


View Commands 


mi_showgrid Show the grid points. 
mi_hidegrid Hide the grid points points. 


mi_grid_snap(’flag’) Setting flag to ’on’ turns on snap to grid, setting flag to ’off’ turns off snap to 
grid. 


mi_setgrid(density,’type’) Change the grid spacing. The density parameter specifies the space 
between grid points, and the type parameter is set to ‘cart’ for cartesian coordinates or ’polar’ for 
polar coordinates. 


mi_refreshview Redraws the current view. 

mi_minimize minimizes the active magnetics input view. 

mi_maximize maximizes the active magnetics input view. 

mi_restore restores the active magnetics input view from a minimized or maximized state. 


mi_resize(width,height) resizes the active magnetics input window client area to width x height. 


Object Properties 


mi_getmaterial(’materialname’) fetches the material specified by materialname from the 
materials library. 


mi_addmaterial(’matname’, mu x, mu y, H c, J, Cduct, Lam d, Phi hmax, lam fill, LamType, Phi 
hx, Phi hy, nstr, dwire) adds a new material with called ’matname’ with the material properties: 


mu x Relative permeability in the x- or r-direction. 


mu y Relative permeability in the y- or z-direction. 


H c Permanent magnet coercivity in Amps/Meter. 


J Applied source current density in Amps/mm2. 


Cduct Electrical conductivity of the material in MS/m. 


Lam d Lamination thickness in millimeters. 


Phi hmax Hysteresis lag angle in degrees, used for nonlinear BH curves. 


— Lam fill Fraction of the volume occupied per lamination that is actually filled with iron (Note 
that this parameter defaults to 1 in the femm preprocessor dialog box because, by default, 
iron completely fills the volume) 


— Lamtype Set to 


* 0 — Not laminated or laminated in plane 
* 1 — laminated x or r 

* 2 — laminated y or z 

x 3 — magnet wire 

* 4 — plain stranded wire 

* 5 — Litz wire 

* 6 — square wire 


Phi hx Hysteresis lag in degrees in the x-direction for linear problems. 


Phi hy Hysteresis lag in degrees in the y-direction for linear problems. 


nstr Number of strands in the wire build. Should be 1 for Magnet or Square wire. 


— dwire Diameter of each of the wire’s constituent strand in millimeters. 


Note that not all properties need be defined — properties that aren’t defined are assigned default 
values. 


mi_addbhpoint(’blockname’,b,h) Adds a B-H data point the the material specified by the 
string ’blockname’. The point to be added has a flux density of b in units of Teslas and a field 
intensity of h in units of Amps/Meter. 


mi_clearbhpoints(’blockname’) Clears all B-H data points associatied with the material 
specified by ’blockname’. 


mi_addpointprop(’pointpropname’,a,j) adds a new point property of name ’pointpropname’ 
with either a specified potential a in units Webers/Meter or a point current j in units of Amps. 
Set the unused parameter pairs to 0. 


mi_addboundprop(’propname’, AO, A1, A2, Phi, Mu, Sig, cO, c1, BdryFormat, ia, oa) 
adds a new boundary property with name ’propname’ 


For a “Prescribed A” type boundary condition, set the AO, A1, A2 and Phi parameters as 
required. Set all other parameters to zero. 


— Fora “Small Skin Depth” type boundary condtion, set the Mu to the desired relative 
permeability and Sig to the desired conductivity in MS/m. Set BdryFormat to 1 and all 
other parameters to zero. 
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To obtain a “Mixed” type boundary condition, set C1 and CO as required and 
BdryFormat to 2. Set all other parameters to zero. 


— Fora “Strategic dual image” boundary, set BdryFormat to 3 and set all other parameters 
to zero. 


— For a “Periodic” boundary condition, set BdryFormat to 4 and set all other parameters 
to zero. 


— Foran “Anti-Perodic” boundary condition, set BdryFormat to 5 set all other parameters 
to zero. 


— Fora “Periodic Air Gap”, set BdryFormat to 6. Parameters ia and oa specify the inner 
and outer boundary angles, respectively. 


— Foran “Anti-periodic Air Gap”, set BdryFormat to 7. The same ia and oa parameters also 


apply here.” 


mi_addcircprop(’circuitname’, i, circuittype) adds a new circuit property with name 
‘circuitname’ with a prescribed current. The circuittype parameter is O for a parallel- 
connected circuit and 1 for a series-connected circuit. 


mi_deletematerial(’materialname’) deletes the material named ’materialname’. 
mi_deleteboundprop(’propname’) deletes the boundary property named ’propname’. 
mi_deletecircuit(’circuitname’) deletes the circuit named circuitname. 
mi_deletepointprop(’pointpropname’) deletes the point property named ’pointpropname’ 


mi_modifymaterial(’BlockName’,propnum,value) This function allows for modification of a 
material’s properties without redefining the entire material (e.g. so that current can be 
modified from run to run). The material to be modified is specified by ’BlockName’. The next 
parameter is the number of the property to be set. The last number is the value to be applied 
to the specified property. The various properties that can be modified are listed below: 








propnum Symbol Description 

0 BlockName Name of the material 

1 Lx x (or r) direction relative permeability 

2 Uy y (or z) direction relative permeability 

3 He Coercivity, Amps/Meter 

4 J Source current density, MA/m? 

5 foj Electrical conductivity, MS/m 

6 diam Lamination thickness, mm 

7 @hmax Hysteresis lag angle for nonlinear problems, degrees 
8 LamFill Iron fill fraction 
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9 LamType 
10 Phx 

degrees 
11 Phy 


0 = None/In plane, 1 = parallel to x, 2=parallel to y 
Hysteresis lag in x-direction for linear problems, 


Hysteresis lag in y-direction for linear problems, 
degrees 





e mi_modifyboundprop(’BdryName’,propnum, value) This function allows for modification of a 
boundary property. The BC to be modified is specified by ’BdryName’. The next parameter is 
the number of the property to be set. The last number is the value to be applied to the 
specified property. The various properties that can be modified are listed below: 





propnum Symbol 


Description 





BdryName 
Ao 
Ai 
A2 
b 


u 
(o 


co 
c1 


O CON DU BP WN FP OO 


BdryFormat 


10 ia 


11 oa 


Name of boundary property 
Prescribed A parameter 

Prescribed A parameter 

Prescribed A parameter 

Prescribed A phase 

Small skin depth relative permeability 
Small skin depth conductivity, MS/m 
Mixed BC parameter 

Mixed BC parameter 

Type of boundary condition: 

0 = Prescribed A 

1 = Small skin depth 


2 = Mixed 
3 = Strategic Dual Image 
4 = Periodic 


5 = Antiperiodic 

6 = Periodic Air Gap 

7 = Antiperiodic Air Gap 

Inner boundary angle for air gap 
element 

Outer boundary angle for air gap 
element 





e mi_modifypointprop(’PointName’,propnum,value) This function allows for modification of a 
point property. The point property to be modified is specified by ’PointName’. The next 
parameter is the number of the property to be set. The last number is the value to be applied 
to the specified property. The various properties that can be modified are listed below: 





propnum Symbol 


Description 
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0 PointName Name of the point property 


1 A Nodal potential, 
Weber/Meter 
2 J Nodal current, Amps 





mi_modifycircprop(’CircName’,propnum,value) This function allows for modification of a 
circuit property. The circuit property to be modified is specified by ’CircName’. The next 
parameter is the number of the property to be set. The last number is the value to be applied 


to the specified property. The various properties that can be modified are listed below: 








propnum Symbol Description 

0 CircName Name of the circuit 
property 

1 i Total current 

2 CircType 0 = Parallel, 1 = Series 





mi_setcurrent(’CircName’,i) sets the current in the circuit specified by ’CircName’ to i. 


4.10 Miscellaneous 


mi_savebitmap( filename’) saves a bitmapped screenshot of the current view to the file 
specified by ’filename’. 


mi_savemetafile(’filename’) saves a metafile screenshot of the current view to the file 
specified by ’filename’. 


mi_refreshview Redraws the current view. 


mi_close Closes current magnetics preprocessor document and destroys magnetics 
preprocessor window. 


mi_shownames(flag) This function allows the user to display or hide the block label 
names on screen. To hide the block label names, flag should be 0. To display the names, 
the parameter should be set to 1. 


mi_readdxf(’filename’) This function imports a dxf file specified by ’filename’. 


mi_savedxf(’filename’) This function saves geometry information in a dxf file specified by 
filename’. 


mi_defineouterspace(Zo,Ro,Ri) defines an axisymmetric external region to be used in 
conjunction with the Kelvin Transformation method of modeling unbounded problems. 
The Zo parameter is the z-location of the origin of the outer region, the Ro parameter is 
the radius of the outer region, and the Ri parameter is the radius of the inner region (i.e. 
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the region of interest). In the exterior region, the permeability varies as a function of 
distance from the origin of the external region. These parameters are necessary to define 
the permeability variation in the external region. 


e mi_attachouterspace marks all selected block labels as members of the external region 
used for modeling unbounded axisymmetric problems via the Kelvin Transformation. 


e mi_detachouterspace undefines all selected block labels as members of the external 
region used for modeling unbounded axisymmetric problems via the Kelvin 
Transformation. 


e mi_attachdefault() marks the selected block label as the default block label. This block 
label is applied to any region that has not been explicitly labeled. 


e mi_detachdefault() undefines the default attribute for the selected block labels. 


e mi_makeABC(n,R,x,y,bc) creates a series of circular shells that emulate the impedance of 
an unbounded domain (i.e. an Improvised Asymptotic Boundary Condition). The n 
parameter contains the number of shells to be used (should be between 1 and 10), R is 
the radius of the solution domain, and (x,y) denotes the center of the solution domain. 
The bc parameter should be specified as 0 for a Dirichlet outer edge or 1 for a Neumann 
outer edge. If the function is called without all the parameters, the function makes up 
reasonable values for the missing parameters. 


e mi_setprevious(filename,prevtype) defines the previous solution to be used as the basis 
for an AC incremental permeability or frozen permeability solution. The prevtype field is 
an integer that specifies whether the solution is to be incremental permeability (1) or 
frozen permeability (2). The filename should include the .ans extension, e.g. 
mi_setprevious(’mymodel.ans’,1) 


5 Magnetics Post Processor Command Set 


There are several scripting commands designed to operate in the postprocessing environment. 


5.1 Data Extraction Commands 


e mo_lineintegral(type) Calculate the line integral for the defined contour 








type name values 1 values 2 values 3 values 4 
0 B.n total B.n avg B.n - - 

1 H.t total H.t avg H.t - - 

2 Contour length surface area - - 
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3 Stress Tensor Force DCr/xforce DCy/zforce 2xr/xforce 2x y/z force 


4 Stress Tensor Torque DC torque 2x torque - - 


5 (B.n)*2 


total (B.n)^2 avg (B.n)*2 - - 





Returns typically two values. The first value is the result of the integral calculation, 
and the second value is the average of the quantity of interest over the contour. 
The only exception is integral 3, which evaluates Maxwell’s stress tensor. This 
integral can return up to four results. For force and torque results, the 2x results 
are only relevant for problems where w!=0. 


mo_blockintegral(type) Calculate a block integral for the selected blocks. The integral types 
are defined as follows: 








Type Definition 

0 AJ 

1 A 

2 Magnetic field energy 

3 Hysteresis and/or lamination losses 

4 Resistive losses 

5 Block cross-section area 

6 Total losses 

7 Total current 

8 Integral of Bx (or B+) over block 

9 Integral of By (or Bz) over block 

10 Block volume 

11 x (or r) part of steady-state Lorentz force 

12 y (or z) part of steady-state Lorentz force 

13 x (or r) part of 2x Lorentz force 

14 y (or z) part of 2x Lorentz force 

15 Steady-state Lorentz torque 

16 2x component of Lorentz torque 

17 Magnetic field coenergy 

18 x (or r) part of steady-state weighted stress tensor 
force 

19 y (or z) part of steady-state weighted stress tensor 
force 

20 x (or r) part of 2x weighted stress tensor force 

21 y (or z) part of 2x weighted stress tensor force 
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22 
23 
24 
25 
26 
27 
28 
29 
30 


Steady-state weighted stress tensor torque 

2x component of weighted stress tensor torque 
R? (i.e. moment of inertia / density) 

x (or r) part of 1x weighted stress tensor force 
y (or z) part of 1x weighted stress tensor force 
1x component of weighted stress tensor torque 
x (or r) part of 1x Lorentz force 

y (or z) part of 1x Lorentz force 

1x component of Lorentz torque 





e mo_getpointvalues(x,y) Get the values associated with the point at (x,y). The function 
returns an array whose contents are, in order: 





Symbol Definition 





A Potential A or flux ọ 

B1 By if planar, Brif axisymmetric 

B2 By if planar, Bzif axisymmetric 

Sig conductivity o 

E stored energy density 

H1 Hy if planar, H-if axisymmetric 

H2 Hy if planar, Hz if axisymmetric 

Je eddy current density 

Js source current density 

Mul Lx if planar, pif axisymmetric 

Mu2 Ly if planar, uz if axisymmetric 

Pe Power density dissipated through ohmic 
losses 

Ph Power density dissipated by hysteresis 

ff Winding fill factor 





The following series of functions retrieves smaller subsets of these results. 


e mo_geta(x,y) Get the potential associated with the point at (x,y). For planar problems, the 
reported potential is vector potential A. For axisymmetric problems, 2mrA is reported. 


e mo_getb(x,y) Get the magnetic flux density associated with the point at (x,y). The return 
value is a list with two elements representing B,and By for planar problems and B,and B; for 
axisymmetric problems. 


e mo_getconductivity(x,y) Gets the conductivity associated with the point at (x,y). 
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mo_getenergydensity(x,y) Gets the magnetic field energy density associated with the point 
at (x,y). 


mo_geth(x,y) Get the magnetic field intensity associated with the point at (x,y). The return 
value is a list with two elements representing Hy and Hy, for planar problems and H; and Hz 
for axisymmetric problems. 


mo_getj(x,y) Get the electric current density associated with the point at (x,y). 


mo_getmu(x,y) Get the relative magnetic permeability associated with the point at (x,y). 
The return value is a list with two elements representing x. and py for planar problems and 
urand uzfor axisymmetric problems. 


mo_getpe(x,y) Get the ohmic loss density associated with the point at (x,y). 


mo_getph(x,y) Get the hysteresis/laminated eddy current loss density associated with the 
point at (x,y). 


mo_getfill(x,y) Get the winding factor (i.e. the average fraction of the volume filled with 
conductor) associated with the point at (x,y). 


mo_makeplot(PlotType,NumPoints,Filename,FileFormat) Allows Octave access to FEMM’s X-Y 
plot routines. If only PlotType or only PlotType and NumPoints are specified, the command is 
interpreted as a request to plot the requested plot type to the screen. If, in addition, the Filename 
parameter is specified, the plot is instead written to disk to the specified file name as an extended 
metafile. If the FileFormat parameter is also, the command is instead interpreted as a command 
to write the data to disk to the specfied file name, rather than display it to make a graphical plot. 
Valid entries for PlotType are: 





PlotType Definition 





Potential 
|B| 

B-n 

B-t 

|H| 

H-n 

H-t 
Jeddy 


CON ONA WN FP OO 


Jsource +Jeddy 





Valid file formats are 





FileFormat Definition 
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Multi-column text with legend 
Multi-column text with no 
legend 

2 Mathematica-style formatting 





For example, if one wanted to plot B-n to the screen with 200 points evaluated to make the graph, 
the command would be: mo_makeplot(2,200) 


If this plot were to be written to disk as a metafile, the command would be: 
mo_makeplot(2,200,’c:\\temp\\myfile.emf’) 
To write data instead of a plot to disk, the command would be of the form: 


mo_makeplot(2,200,’c:\\temp\\myfile.txt’,O) 


mo_getprobleminfo Returns info on problem description. Returns four values: 





Value Definition 





1 problem type 

2 frequency in Hz 

3 problem depth in meters 

4 length unit used to draw the problem, represented in meters 





mo_getcircuitproperties(’circuit’) Used primarily to obtain impedance information associated 
with circuit properties. Properties are returned for the circuit property named ’circuit’. Six values 
are returned by the function. In order, these parameters are: 

— current Current carried by the circuit. 


— volts Voltage drop across the circuit in the circuit. 
— flux Circuit’s flux linkage 
mo_getgapb("BdryName", angle) Computes the radial and tangential flux density on the centerline 


of the specified air gap element name (BdryName) at the specified angle. The angle is specified in 
degrees. 


mo_getgapa("BdryName",angle) Computes the magnetic vector potential on the centerline of the 
specified air gap element name (BdryName) at the specified angle. The angle is specified in 
degrees. 


mo_gapintegral("BdryName",inttype) Computes an integral specifiedy by inttype over an air gap 
element specified by BdryName. Values for inttype in are: 





0 DC Torque 
1 DC Force 
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5.2 


Stored Energy 

2X Torque 

2X Force 
Interaction Torque 
Interaction Force 


nu AUN 





Torque and energy integrals return one result; force integrals return two results for the x- and y- 
direction force, respectively. 


mo_getgapharmonics("BdryName",n) Returns the acc, acs, brc, brs, btc, and bts. These quantities 
represent the components of vector potential, radial flux density, and tangential flux density on 
the centerline of the specified air gap element at the specified angle. For the n harmonic, vector 
potential, radial flux density, and tangential flux density can be represented explicitly as functions 
of angle via: 


A acc*cos(n@)+acs*sin(n8) 
B, = brc*cos(n8)+brs*sin(n8) 
Bə = btc*cos(n8)+bts*sin(n8) 


If the function is called with just the BdryName, the function returns the number of harmonics 
available. 


Selection Commands 


mo_seteditmode(mode) Sets the mode of the postprocessor to point, contour, or area mode. Valid 
entries for mode are ’point’, ‘contour’, and ‘area’. 


mo_selectblock(x,y) Select the block that contains point (x,y). 


mo_groupselectblock(n) Selects all the blocks that are labeled by block labels that are members of 
group n. If no number is specified (i.e. mo_groupselectblock() ), all blocks are selected. 


mo_addcontour(x,y) Adds a contour point at (x,y). If this is the first point then it starts a contour, 
if there are existing points the contour runs from the previous point to this point. The 
mo_addcontour command has the same functionality as a right-button-click contour point 
addition when the program is running in interactive mode. 


mo_bendcontour(angle,anglestep) Replaces the straight line formed by the last two points in the 
contour by an arc that spans angle degrees. The arc is composed of many straight lines, each of 
which is constrained to span no more than anglestep degrees. The angle parameter can take on 
values from -180 to 180 degrees. The anglestep parameter must be greater than zero. If there are 
less than two points defined in the contour, this command is ignored. 
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e mo_selectpoint(x,y) Adds a contour point at the closest input point to (x,y). If the selected point 
and a previously selected point lie at the ends of an arcsegment, a contour is added that traces 
along the arcsegment. The mo_selectpoint command has the same functionality as the left- 
button-click contour point selection when the program is running in interactive mode. 


e mo_clearcontour Clear a prevously defined contour 


e mo_clearblock Clear block selection 


5.3. Zoom Commands 


e mo_zoomnatural Zoom to the natural boundaries of the geometry. 
e mo_zoomin Zoom in one level. 
e mo_zoomout Zoom out one level. 


e mo_zoom(x1,y1,x2,y2) Zoom to the window defined by lower left corner (x1,y1) and upper right 
corner (x2,y2). 


5.4 View Commands 


mo_showmesh Show the mesh. 

mo_hidemesh Hide the mesh. 

mo_showpoints Show the node points from the input geometry. 
mo_hidepoints Hide the node points from the input geometry. 


mosmooth(’flag’) This function controls whether or not smoothing is applied to the B and H fields, 
which are naturally piece-wise constant over each element. Setting flag equal to ’on’ turns on 
smoothing and setting flag to ‘off’ turns off smoothing. 


mo_showgrid Show the grid points. 
mo_hidegrid Hide the grid points points. 


mo_grid_snap(’flag’) Setting flag to ’on’ turns on snap to grid, setting flag to ’off’ turns off snap to 
grid. 


mo_setgrid(density,’type’) Change the grid spacing. The density parameter specifies the space 
between grid points, and the type parameter is set to cart’ for cartesian coordinates or ’polar’ for 
polar coordinates. 


mo_hidedensityplot hides the flux density plot. 
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e mo_showdensityplot(legend,gscale,upper_B,lower_B,type) Shows the flux density plot with options: 


— legend Set to 0 to hide the plot legend or 1 to show the plot legend. 

—  gscale Set to 0 for a colour density plot or 1 for a grey scale density plot. 

— _upper_B Sets the upper display limit for the density plot. — lower_B Sets the lower display limit 
for the density plot. 


— type Type of density plot to display. Valid entries are ’mag’, ’real’, and ‘imag’ for magnitude, real 
component, and imaginary component of B, respectively. Alternatively, current density can be 
displayed by specifying ’jmag’, ’jreal’, and ’jimag’ for magnitude, real component, and imaginary 
component of J, respectively. 


e mo_hidecontourplot Hides the contour plot. 


e mo_showcontourplot(numcontours,lower_A,upper_A,type) shows the A contour plot with options: 


— numcontours Number of A equipotential lines to be plotted. 
— upper_A Upper limit for A contours. 
— lower_A Lower limit for A contours. 


— type Choice of ’real’, ’imag’, or ’both’ to show either the real, imaginary of both real and 
imaginary components of A. 


mo_showvectorplot(type,scalefactor) controls the display of vectors denoting the field strength and 
direction. The parameters taken are the type of plot, which should be set to O for no vector plot, 1 
for the real part of flux density B; 2 for the real part of field intensity H; 3 for the imaginary part of B; 
4 for the imaginary part of H; 5 for both the real and imaginary parts of B; and 6 for both the real and 
imaginary parts of H. The scalefactor determines the relative length of the vectors. If the scale is set 
to 1, the length of the vectors is chosen so that the highest flux density corresponds to a vector that 
is the same length as the current grid size setting. 


e mo_minimize minimizes the active magnetics output view. 
e mo_maximize maximizes the active magnetics output view. 
e mo_restore restores the active magnetics output view from a minimized or maximized state. 


e mo_resize(width,height) resizes the active magnetics output window client area to width x height. 


5.5 Miscellaneous 


e mo_close Closes the current post-processor instance. 


e mo_refreshview Redraws the current view. 
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mo_reload Reloads the solution from disk. 


mo_savebitmap(’filename’) saves a bitmapped screen shot of the current view to the file specified 
by ’filename’. 


mo_savemetafile(’filename’) saves a metafile screenshot of the current view to the file specified 
by ’filename’. 


mo_shownames(flag) This function allows the user to display or hide the block label names on 
screen. To hide the block label names, flag should be 0. To display the names, the parameter should 
be set to 1. 


mo_numnodes Returns the number of nodes in the in-focus magnetics output mesh. 


mo_numelements Returns the number of elements in the in-focus magnets output mesh. 


mo_getnode(n) Returns the (x,y) or (r,z) position of the nth mesh node. 


mo_getelement(n) returns the following proprerties for the nth element: 


. Index of first element node 

. Index of second element node 

. Index of third element node 

. x (or r) coordinate of the element centroid 
. y (or z) coordinate of the element centroid 


. element area using the length unit defined for the problem 


N DOD WU BP WN e 


. group number associated with the element 


Electrostatics Preprocessor Command Set 


Several different commands are available in the preprocessor. Two naming conventions can be used: 
one which separates words in the command names by underscores, and one that eliminates the 
underscores. 


6.1 


Object Add/Remove Commands 


ei_addnode(x,y) Add a new node at x,y 


ei_addsegment(x1,y1,x2,y2) Add a new line segment from node closest to (x1,y1) to node closest 
to (x2,y2) 


ei_addblocklabel(x,y) Add a new block label at (x,y) 
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6.2 


ei_addarc(x1,y1,x2,y2,angle,maxseg) Add a new arc segment from the nearest node to (x1,y1) to 
the nearest node to (x2,y2) with angle ‘angle’ divided into ‘maxseg’ segments. 


ei_drawline(x1,y1,x2,y2) Adds nodes at (x1,y1) and (x2,y2) and adds a line between the nodes. 


ei_drawarc(x1,y1,x2,y2,angle,maxseg) Adds nodes at (x1,y1) and (x2,y2) and adds an arc of the 
specified angle and discretization connecting the nodes. 


ei_drawrectangle(x1,y1,x2,y2) Adds nodes at the corners of a rectangle defined by the points 
(x1,y1) and (x2,y2), then adds segments connecting the corners of the rectangle. 


ei_deleteselected Delete all selected objects. 
ei_deleteselectednodes Delete selected nodes. 
ei_deleteselectedlabels Delete selected block labels. 
ei_deleteselectedsegments Delete selected segments. 


ei_deleteselectedarcsegments Delete selects arcs. 


Geometry Selection Commands 


ei_clearselected Clear all selected nodes, blocks, segments and arc segments. 
ei_selectsegment(x,y) Select the line segment closest to (x,y) 

ei_selectnode(x,y) Select the node closest to (x,y). Returns the coordinates of the selected node. 
ei_selectlabel(x,y) Select the label closet to (x,y). Returns the coordinates of the selected label. 
ei_selectarcsegment(x,y) Select the arc segment closest to (x,y) 


ei_selectgroup(n) Select the n* group of nodes, segments, arc segments and block labels. 
This function will clear all previously selected elements and leave the edit mode in 4 (group) 


ei_selectcircle(x,y,R,editmode) selects objects within a circle of radius R centered at (x,y). If only x, 
y, and R parameters are given, the current edit mode is used. If the editmode parameter is used, 
0 denotes nodes, 2 denotes block labels, 2 denotes segments, 3 denotes arcs, and 4 specifies that 
all entity types are to be selected. 


ei_selectrectangle(x1,y1,x2,y2,editmode) selects objects within a rectangle defined by points 
(x1,y1) and (x2,y2). If no editmode parameter is supplied, the current edit mode is used. If the 
editmode parameter is used, 0 denotes nodes, 2 denotes block labels, 2 denotes segments, 3 
denotes arcs, and 4 specifies that all entity types are to be selected. 
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6.3 


Object Labeling Commands 


e ej_setnodeprop(’propname’,groupno, ’inconductor’) Set the selected nodes to have the nodal 


property ’propname’ and group number groupno. The ’inconductor’ string specifies which 
conductor the node belongs to. If the node doesn’t belong to a named conductor, this parameter 
can be set to ’<None>’. 


ei_setblockprop(’blockname’, automesh, meshsize, group) Set the selected block labels to have 


the properties: Block property ’blockname’. 


automesh: O = mesher defers to mesh size constraint defined in meshsize, 1 = mesher 
automatically chooses the mesh density. meshsize: size constraint on the mesh in the block 


marked by this label. A member of group number group 

ei_setsegmentprop(’name’, elmsize, automesh, hide, group, ’inconductor’) Set the select 
segments to have: 

Boundary property ’name’ 

Local element size along segment no greater than elmsize 


automesh: 0 = mesher defers to the element constraint defined by elementsize, 1 = mesher 
automatically chooses mesh size along the selected segments hide: O = not hidden in post- 
processor, 1 == hidden in post processor 


A member of group number group 


A member of the conductor specified by the string ’inconductor’. If the segment is not part of a 
conductor, this parameter can be specified as ‘<None>’. 


ei_setarcsegmentprop(maxsegdeg, ’propname’, hide, group, ’inconductor’) Set the selected arc 
segments to: 


Meshed with elements that span at most maxsegdeg degrees per element 


Boundary property ’propname’ 
hide: 0 = not hidden in post-processor, 1 == hidden in post processor 
A member of group number group 


A member of the conductor specified by the string ’inconductor’. If the segment is not part of a 
conductor, this parameter can be specified as ’*<None>’. 


e ei_setgroup(n) Set the group associated of the selected items to n. 
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6.4 


Problem Commands 


e ei_probdef(units,type,precision,depth,minangle) changes the problem definition. The units 


6.5 


6.6 


parameter specifies the units used for measuring length in the problem domain. 

Valid ’units’ entries are ‘inches’, ’millimeters’, ’centimeters’, mils’, ’meters, and ’micrometers’. 
Set problemtype to ’planar’ for a 2-D planar problem, or to ’axi’ for an axisymmetric problem. 
The precision parameter dictates the precision required by the solver. For example, entering 1.E- 
8 requires the RMS of the residual to be less than 10-8. The fourth parameter represents the 
depth of the problem in the into-the-page direction for 2-D planar problems-—just enter 0 as a 
placeholder for axisymmetric problems. The sixth parameter represents the minimum angle 
constraint sent to the mesh generator (usually 30 degrees). 


ei_analyze(flag) runs the electrostatics solver. The flag parameter controls whether the solver 
window is visible or minimized. For a visible window, specify 0. For a minimized window, flag 
should be set to 1. If no value is specified for flag, the visibility of the solver is inherited from the 
main window, i.e. if the main window is minimized, the solver runs minimized, too. 


ei_loadsolution loads and displays the solution corresponding to the current geometry. 


ei_setfocus(’documentname’) Switches the electrostatics input file upon which commands are to 
act. If more than one electrostatics input file is being edited at a time, this command can be used 
to switch between files so that the mutiple files can be operated upon programmatically. 
’documentname’ should contain the name of the desired document as it appears on the window’s 
title bar. 


ei_saveas(’filename’) saves the file with name ’filename’. 


Mesh Commands 


ei_createmesh runs triangle to create a mesh. Note that this is not a necessary precursor of 
performing an analysis, as ei_analyze will make sure the mesh is up to date before running an 
analysis. The number of elements in the mesh is pushed back onto the lua stack. 


ei_showmesh toggles the flag that shows or hides the mesh. 


ei_purgemesh clears the mesh out of both the screen and memory. 


Editing Commands 


ei_copyrotate(bx, by, angle, copies ) 


— bx, by — base point for rotation 


— angle — angle by which the selected objects are incrementally shifted to make each copy. 
angle is measured in degrees. 
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— copies — number of copies to be produced from the selected objects. 
e ei_copyrotate2(bx, by, angle, copies, editaction ) 


— bx, by — base point for rotation 


angle — angle by which the selected objects are incrementally shifted to make each copy. 
angle is measured in degrees. 


— copies — number of copies to be produced from the selected objects. 


editaction 0 —nodes, 1 — lines (segments), 2 —block labels, 3 — arc segments, 4- group 
e ei_copytranslate(dx, dy, copies) 


— dx,dy — distance by which the selected objects are incrementally shifted. 


— copies — number of copies to be produced from the selected objects. 
e ei_copytranslate2(dx, dy, copies, editaction) 


— dx,dy — distance by which the selected objects are incrementally shifted. 
— copies — number of copies to be produced from the selected objects. 


— editaction 0 —-nodes, 1 — lines (segments), 2 -block labels, 3 — arc segments, 4- group 
e ei_createradius(x,y,r) turns a corner located at (x,y) into a curve of radius r. 
e eji_moverotate(bx,by,shiftangle) 


— bx, by — base point for rotation 


— shiftangle — angle in degrees by which the selected objects are rotated. 
e ei_moverotate2(bx,by,shiftangle, editaction) 


— bx, by — base point for rotation 
— shiftangle — angle in degrees by which the selected objects are rotated. 


— editaction 0 —-nodes, 1 — lines (segments), 2 -block labels, 3 — arc segments, 4- group 
e ei_movetranslate(dx,dy) 

— dx,dy — distance by which the selected objects are shifted. 
e eji_movetranslate2(dx,dy,editaction) 


— dx,dy — distance by which the selected objects are shifted. 


— editaction 0 —-nodes, 1 — lines (segments), 2 -block labels, 3 — arc segments, 4- group 
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6.7 


6.8 


ei_scale(bx, by,scalefactor) 


— bx, by — base point for scaling 


— scalefactor — a multiplier that determines how much the selected objects are scaled 
ei_scale2(bx,by,scalefactor,editaction) 


— bx, by — base point for scaling 
— scalefactor — a multiplier that determines how much the selected objects are scaled 


— editaction 0 —nodes, 1 — lines (segments), 2 —block labels, 3 — arc segments, 4- group 


ei_mirror(x1,y1,x2,y2) mirror the selected objects about a line passing through the points (x1,y1) 
and (x2,y2). 


ei_mirror2(x1,y1,x2,y2,editaction) mirror the selected objects about a line passing through the 
points (x1,y1) and (x2,y2). Valid editaction entries are O for nodes, 1 for lines (segments), 2 for 
block labels, 3 for arc segments, and 4 for groups. 


ei_seteditmode(editmode) Sets the current editmode to: 
— ’nodes’ - nodes 
— ’segments’ - line segments 
— ’arcsegments’ - arc segments 
— ’blocks’ - block labels 


— ’group’ - selected group 


This command will affect all subsequent uses of the other editing commands, if they are used 
without the editaction parameter. 


Zoom Commands 


|” 


ei_zoomnatural zooms to a “natural” view with sensible extents. 
ei_zoomout zooms out by a factor of 50%. 
ei_zoomin zoom in by a factor of 200%. 


ei_zoom(x1,y1,x2,y2) Set the display area to be from the bottom left corner specified by (x1,y1) to 
the top right corner specified by (x2,y2). 


View Commands 


ei_showgrid Show the grid points. 
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6.9 


ei_hidegrid Hide the grid points points. 


ei_gridsnap(’flag’) Setting flag to “on” turns on snap to grid, setting flag to “off” turns off snap to 
grid. 


ei_setgrid(density,’type’) Change the grid spacing. The density parameter specifies the space 
between grid points, and the type parameter is set to ‘cart’ for Cartesian coordinates or ’polar’ for 
polar coordinates. 


ei_refreshview Redraws the current view. 

ei_minimize minimizes the active electrostatics input view. 

ei_maximize maximizes the active electrostatics input view. 

ei_restore restores the active electrostatics input view from a minimized or maximized state. 


ei_resize(width,height) resizes the active electrostatics input window client area to width x height. 


Object Properties 
ei_getmaterial(’materialname’) fetches the material specified by materialname from the materials 
library. 
ei_addmaterial(’matname’, ex, ey, qv) adds a new material with called ’matname’ with the 
material properties: 


ex Relative permittivity in the x- or r-direction. 
ey Relative permittivity in the y- or z-direction. 
qv Volume charge density in units of C/m? 


ei_addpointprop(’pointname’,Vp,qp) adds a new point property of name ’pointname’ with either 
a specified potential Vp a point charge density qp in units of C/m. . 


ei_addboundprop(’boundname’, Vs, qs, cO, c1, BdryFormat) adds a new boundary property with 
name ’boundname’ 


For a “Fixed Voltage” type boundary condition, set the Vs parameter to the desired voltage and all 
other parameters to zero. 


To obtain a “Mixed” type boundary condition, set C1 and CO as required and BdryFormat to 
1. Set all other parameters to zero. 


To obtain a prescribes surface charge density, set qs to the desired charge density in C/m? and set 
BdryFormat to 2. 
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For a “Periodic” boundary condition, set BdryFormat to 3 and set all other parameters to zero. 
For an “Anti-Perodic” boundary condition, set BdryFormat to 4 set all other parameters to zero. 
ei_addconductorprop(’conductorname’, Vc, qc, conductortype) adds a new conductor property 
with name ’conductorname’ with either a prescribed voltage or a prescribed total charge. Set the 


unused property to zero. The conductortype parameter is O for prescribed charge and 1 for 
prescribed voltage. 


ei_deletematerial(’materialname’) deletes the material named ’materialname’. 
ei_deleteboundprop(’boundname’) deletes the boundary property named ’boundname’. 
ei_deleteconductor(’conductorname’) deletes the conductor named conductorname. 
ei_deletepointprop(’pointname’) deletes the point property named ’pointname’ 


ei_modifymaterial(’BlockName’,propnum,value) This function allows for modification of a 
material’s properties without redefining the entire material (e.g. so that current can be modified 
from run to run). The material to be modified is specified by ’BlockName’. The next parameter is 
the number of the property to be set. The last number is the value to be applied to the specified 
property. The various properties that can be modified are listed below: 





propnum Symbol Description 

0 BlockName Name of the material 

1 ex x (or r) direction relative permittivity 
2 ey y (or z) direction relative permittivity 
3 qs Volume charge 


ei_modifyboundprop(’BdryName’,propnum,value) This function allows for modification of a 
boundary property. The BC to be modified is specified by ’BdryName’. The next parameter is the 
number of the property to be set. The last number is the value to be applied to the specified 
property. The various properties that can be modified are listed below: 





propnum Symbol Description 

0 BdryName Name of boundary property 
1 Vs Fixed Voltage 

2 qs Prescribed charge density 

3 c0 Mixed BC parameter 

4 c1 Mixed BC parameter 
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5 BdryFormat Type of boundary condition: 
O = Prescribed V 


1 = Mixed 

2 = Surface charge density 
3 = Periodic 

4 = Antiperiodic 


ei_modifypointprop(’PointName’,propnum,value) This function allows for modification of a point 
property. The point property to be modified is specified by ’PointName’. The next parameter is the 
number of the property to be set. The last number is the value to be applied to the specified 
property. The various properties that can be modified are listed below: 





propnum Symbol Description 
0 PointName Name of the point property 
1 Vp Prescribed nodal voltage 

2 qp Point charge density in C/m 


ei_modifyconductorprop(’ConductorName’,propnum,value) This function allows for modification 
of a conductor property. The conductor property to be modified is specified by ’ConductorName’. 
The next parameter is the number of the property to be set. The last number is the value to be 
applied to the specified property. The various properties that can be modified are listed below: 





propnum Symbol Description 

0 ConductorName Name of the conductor property 

1 Vc Conductor voltage 

2 qc Total conductor charge 

3 ConductorType 0 = Prescribed charge, 1 = Prescribed voltage 


6.10 Miscellaneous 


ei_savebitmap( filename’) saves a bitmapped screenshot of the current view to the file specified 
by ’filename’. 


ei_savemetafile(’filename’) saves a metafile screenshot of the current view to the file specified by 
filename’. 


ei_refreshview Redraws the current view. 
ei_close closes the preprocessor window and destroys the current document. 


ei_shownames(flag) This function allows the user to display or hide the block label names on 
screen. To hide the block label names, flag should be 0. To display the names, the parameter should 
be set to 1. 
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ei_readdxf(’filename’) This function imports a dxf file specified by ’filename’. 


ei_savedxf(’filename’) This function saves geometry information in a dxf file specified by 
filename’. 


ei_defineouterspace(Zo,Ro,Ri) defines an axisymmetric external region to be used in conjuction 
with the Kelvin Transformation method of modeling unbounded problems. The Zo parameter is 
the z-location of the origin of the outer region, the Ro parameter is the radius of the outer region, 
and the Ri parameter is the radius of the inner region (i.e. the region of interest). In the exterior 
region, the permeability varies as a function of distance from the origin of the external region. 
These parameters are necessary to define the permeability variation in the external region. 


ei_attachouterspace marks all selected block labels as members of the external region used for 
modeling unbounded axisymmetric problems via the Kelvin Transformation. 


ei_detachouterspace undefines all selected block labels as members of the external region used 
for modeling unbounded axisymmetric problems via the Kelvin Transformation. 


ei_attachdefault() marks the selected block label as the default block label. This block label is 
applied to any region that has not been explicitly labeled. 


ei_detachdefault() undefines the default attribute for the selected block labels. 


ei_makeABC(n,R,x,y,bc) creates a series of circular shells that emulate the impedance of an 
unbounded domain (i.e. an Inprovised Asymptotic Boundary Condition). The n parameter contains 
the number of shells to be used (should be between 1 and 10), R is the radius of the solution 
domain, and (x,y) denotes the center of the solution domain. The bc parameter should be specified 
as 0 for a Dirichlet outer edge or 1 for a Neumann outer edge. If the function is called without all 
the parameters, the function makes up reasonable values for the missing parameters. 


Electrostatics Post Processor Command Set 


There are several scripting commands designed to operate in the postprocessor. As with the 
preprocessor commands, these commands can be used with either the underscore naming or with the 
no-underscore naming convention. 


7.1 


Data Extraction Commands 


e eo_lineintegral(type) Calculates the line integral for the defined contour type Integral 





0 E-t 
1 D-n 
2 Contour length 
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3 Force from stress tensor 
4 Torque from stress tensor 


This integral a single result for field intensity, flux density, contour length, and torque. For force, 
an array with two elements is returned representing force in the x and y or r and z directions. 


e eo_blockintegral(type) Calculates a block integral for the selected blocks 


type Integral 





Stored Energy 

Block Cross-section 

Block Volume 

Average D over the block 4 Average E over the block 
Weighted Stress Tensor Force 


Ua wU N BeO 


6 Weighted Stress Tensor Torque 


This integral a single result for energy, volume and torque. For force, D and E£, the function 
returns an array with two elements is returned representing components in the x and y or rand 
z directions. 


e eo_getpointvalues(x,y) Gets the values associated with the point at (x,y). The function returns 
an array of results whose elements represent: 


Symbol Definition 





V Voltage 

Dx x- or r- direction component of displacement 

Dy y- or z- direction component of displacement 

Ex x- or r- direction component of electric field intensity 
Ey y- or z- direction component of electric field intensity 
ex x- or r- direction component of permittivity 

ey y- or z- direction component of permittivity 

nrg electric field energy density 


e eo_getv(x,y) Gets the voltage associated with the point at (x,y). 


e eo _getd(x,y) Gets the electric flux density associated with the point at (x,y). The return value is 
a list with two elements representing Dx and Dy for planar problems and D; and Dz for 
axisymmetric problems. 


e eo_gete(x,y) Gets the electric field intensity associated with the point at (x,y). The return value 
is a list with two element representing Ex and Ey for planar problems and E, and Ez for 
axisymmetric problems. 
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e eo_getenergydensity(x,y) Gets the electric field energy density associated with the point at (x,y). 


e eo_getperm(x,y) Gets the electric permittivity associated with the point at (x,y). The return value 
is a list with two elements representing €x and £y for planar problems and €r and ez for 
axisymmetric problems. 


e eo_makeplot(PlotType,NumPoints, Filename,FileFormat) Allows Octave to access to the X-Y plot 
routines. If only PlotType or only PlotType and NumPoints are specified, the command is 
interpreted as a request to plot the requested plot type to the screen. If, in addition, the 
Filename parameter is specified, the plot is instead written to disk to the specified file name as 
an extended metafile. If the FileFormat parameter is also, the command is instead interpreted 
as a command to write the data to disk to the specfied file name, rather than display it to make 
a graphical plot. Valid entries for PlotType are: 


PlotType Definition 





V (Voltage) 

|D| (Magnitude of flux density) 
D . n (Normal flux density) 

D . t (Tangential flux density) 

JE] (Magnitude of field intensity) 
E . n (Normal field intensity) 


nu A WN FP OO 


E . t (Tangential field intensity) 


Valid file formats are: 


FileFormat Definition 





Multi-column text with legend 
Multi-column text with no legend 
Mathematica-style formatting 


For example, if one wanted to plot V to the screen with 200 points evaluated to make the graph, 
the command would be: eo_makeplot(0,200) 

If this plot were to be written to disk as a metafile, the command would be: 
eo_makeplot(0,200,’c:\\temp\\temp.emf’) 

To write data instead of a plot to disk, the command would be of the form: 
eo_makeplot(0,200,’c:\\temp\\temp.txt’,0) 


e eo _getprobleminfo Returns info on problem description. Returns three values: the Problem type 
(0 for planar and 1 for axisymmetric); the depth assumed for planar problems in units of meters; 
and the length unit used to draw the problem, represented in meters. 
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7.2 


7.3 


eo_getconductorproperties(’conductor’) Properties are returned for the conductor property 
named ’conductor’. Two values are returned: The voltage of the specified conductor, and the 
charge carried on the specified conductor. 


Selection Commands 


eo_seteditmode(mode) Sets the mode of the postprocessor to point, contour, or area mode. 
Valid entries for mode are ’point’, ’contour’, and ’area’. 


eo_selectblock(x,y) Select the block that contains point (x,y). 


eo_groupselectblock(n) Selects all of the blocks that are labeled by block labels that are members 
of group n. If no number is specified (i.e. eo_groupselectblock), all blocks are selected. 


eo_selectconductor(’name’) Selects all nodes, segments, and arc segments that are part of the 
conductor specified by the string (‘name’). This command is used to select conductors for the 
purposes of the “weighted stress tensor” force and torque integrals, where the conductors are 
points or surfaces, rather than regions (i.e. can’t be selected with eo_selectblock). 


eo_addcontour(x,y) Adds a contour point at (x,y). If this is the first point then it starts a contour, 
if there are existing points the contour runs from the previous point to this point. The 
eo_addcontour command has the same functionality as a right-button-click contour point 
addition when the program is running in interactive mode. 


eo_bendcontour(angle,anglestep) Replaces the straight line formed by the last two points in the 
contour by an arc that spans angle degrees. The arc is composed of many straight lines, each of 
which is constrained to span no more than anglestep degrees. The angle parameter can take on 
values from -180 to 180 degrees. The anglestep parameter must be greater than zero. If there are 
less than two points defined in the contour, this command is ignored. 


eo_selectpoint(x,y) Adds a contour point at the closest input point to (x,y). If the selected point 
and a previously selected point lie at the ends of an arcsegment, a contour is added that traces 
along the arcsegment. The selectpoint command has the same functionality as the left-button- 
click contour point selection when the program is running in interactive mode. 


eo_clearcontour Clear a prevously defined contour 


eo_clearblock Clear block selection 


Zoom Commands 


e0_zoomnatural Zoom to the natural boundaries of the geometry. 


eo_zoomin Zoom in one level. 
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7.4 


eo_zoomout Zoom out one level. 


eo_zoom(x1,y1,x2,y2) Zoom to the window defined by lower left corner (x1,y1) and upper right 
corner (x2,y2). 


View Commands 
eo_showmesh Show the mesh. 
eo_hidemesh Hide the mesh. 
eo_showpoints Show the node points from the input geometry. 
eo_hidepoints Hide the node points from the input geometry. 


eo_smooth(’flag’) This function controls whether smoothing is applied to the D and E fields which 
are naturally piece-wise constant over each element. Setting flag equal to ’on’ turns on smoothing 
and setting flag to ’off’ turns off smoothing. 


eo_showgrid Show the grid points. 


eo_hidegrid Hide the grid points points. 
eo_gridsnap(’flag’) Setting flag to “on” turns on snap to grid, setting flag to ’off’ turns off snap 


to grid. 


eo_setgrid(density,’type’) Change the grid spacing. The density parameter specifies the space 
between grid points, and the type parameter is set to ‘cart’ for Cartesian coordinates or ‘polar’ for 
polar coordinates. 


eo_hidedensityplot hides the flux density plot. 


eo_showdensityplot(legend,gscale,type,upper D,lower D) Shows the flux density plot with options: 
legend Set to 0 to hide the plot legend or 1 to show the plot legend. 

gscale Set to 0 for a colour density plot or 1 for a grey scale density plot. 

upper D Sets the upper display limit for the density plot. lower D Sets the 

lower display limit for the density plot. 

type Sets the type of density plot. A value of 0 plots voltage, 1 plots the magnitude of D, and 2 
plots the magnitude of E 


eo_hidecontourplot Hides the contour plot. 


eo_showcontourplot(numcontours,lower V,upper V) shows the V contour plot with options: 
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7.5 


numcontours Number of equipotential lines to be plotted. 
upper V Upper limit for contours. lower V Lower limit for 


contours. 


eo_showvectorplot(type,scalefactor) controls the display of vectors denoting the field strength 
and direction. The parameters taken are the type of plot, which should be set to O for no vector 
plot, 1 for flux density D, and 2 for field intensity E. The scalefactor determines the relative length 
of the vectors. If the scale is set to 1, the length of the vectors is chosen so that the highest flux 
density corresponds to a vector that is the same length as the current grid size setting. 


eo_minimize minimizes the active electrostatics output view. 

eo_maximize maximizes the active electrostatics output view. 

eo_restore restores the active electrostatics output view from a minimized or maximized state. 
eo_resize(width,height) resizes the active electrostatics output window client area to width x 


height. 


Miscellaneous 


eo_close close the current postprocessor window. 
eo_refreshview Redraws the current view. 
eo_reload Reloads the solution from disk. 


eo_savebitmap(’filename’) saves a bitmapped screen shot of the current view to the file specified 
by ’filename’. 


eo_savemetafile(’filename’) saves a metafile screenshot of the current view to the file specified 
by ’filename’. 


eo_shownames(flag) This function allows the user to display or hide the block label names on 
screen. To hide the block label names, flag should be 0. To display the names, the parameter should 
be set to 1. 


eo_numnodes Returns the number of nodes in the in-focus electrostatics output mesh. 
eo_numelements Returns the number of elements in the in-focus electrostatics output mesh. 
eo_getnode(n) Returns the (x,y) or (r,z) position of the nth mesh node. 


eo_getelement(n) returns the following properties for the nth element: 
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. Index of first element node 

. Index of second element node 

. Index of third element node 

. x (or r) coordinate of the element centroid 
. y (or z) coordinate of the element centroid 


. element area using the length unit defined for the problem 


N DO UW BP WN PR 


. group number associated with the element 


Heat Flow Preprocessor Lua Command Set 


Several different commands are available in the preprocessor. Two naming conventions can be 
used: one which separates words in the command names by underscores, and one that eliminates 
the underscores. 


8.1 


8.2 


Object Add/Remove Commands 
hi_addnode(x,y) Add a new node at x,y 


hi_addsegment(x1,y1,x2,y2) Add a new line segment from node closest to (x1,y1) to node closest 
to (x2,y2) 


hi_addblocklabel(x,y) Add a new block label at (x,y) 


hi_addarc(x1,y1,x2,y2,angle, maxseg) Add a new arc segment from the nearest node to (x1,y1) to 
the nearest node to (x2,y2) with angle ‘angle’ divided into ‘maxseg’ segments. 


hi_deleteselected Delete all selected objects. 
hi_deleteselectednodes Delete selected nodes. 
hi_deleteselectedlabels Delete selected block labels. 
hi_deleteselectedsegments Delete selected segments. 


hi_deleteselectedarcsegments Delete selects arcs. 


Geometry Selection Commands 


hi_clearselected() Clear all selected nodes, blocks, segments and arc segments. 
hi_selectsegment(x,y) Select the line segment closest to (x,y) 


hi_selectnode(x,y) Select the node closest to (x,y). Returns the coordinates of the selected node. 
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8.3 


hi_selectlabel(x,y) Select the label closet to (x,y). Returns the coordinates of the selected label. 
hi_selectarcsegment(x,y) Select the arc segment closest to (x,y) 


hi_selectgroup(n) Select the n“ group of nodes, segments, arc segments and block labels. This 
function will clear all previously selected elements and leave the edit mode in 4 (group) 


hi_selectcircle(x,y,R,editmode) selects objects within a circle of radius R centered at (x,y). If only 
x, y, and R paramters are given, the current edit mode is used. If the editmode parameter is used, 
0 denotes nodes, 2 denotes block labels, 2 denotes segments, 3 denotes arcs, and 4 specifies 
that all entity types are to be selected. 


hi_selectrectangle(x1,y1,x2,y2,editmode) selects objects within a rectangle defined by points 
(x1,y1) and (x2,y2). If no editmode parameter is supplied, the current edit mode is used. If the 
editmode parameter is used, 0 denotes nodes, 2 denotes block labels, 2 denotes segments, 3 
denotes arcs, and 4 specifies that all entity types are to be selected. 


Object Labeling Commands 


hi_setnodeprop("propname",groupno, "inconductor") Set the selected nodes to have 
the nodal property "propname" and group number groupno. The "inconductor" string 
specifies which conductor the node belongs to. If the node doesn’t belong to a named 
conductor, this parameter can be set to "<None>". 


hi_setblockprop("blockname", automesh, meshsize, group) Set the selected block labels to have 


the properties: Block property "blockname". 


automesh: 0 = mesher defers to mesh size constraint defined in meshsize, 1 = mesher 
automatically chooses the mesh density. meshsize: size constraint on the mesh in the 


block marked by this label. A member of group number group 
hi_setsegmentprop("propname", elementsize, automesh, hide, group, "inconductor") Set the 
select segments to have: 

Boundary property "propname" 

Local element size along segment no greater than elementsize 


automesh: O = mesher defers to the element constraint defined by elementsize, 1 = 
mesher automatically chooses mesh size along the selected segments hide: 0 = not 
hidden in post-processor, 1 == hidden in post processor 


A member of group number group 


A member of the conductor specified by the string "inconductor". If the segment is not part 
of a conductor, this parameter can be specified as "<None>". 
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e hi_setarcsegmentprop(maxsegdeg, "propname", hide, group, "inconductor") Set the selected 


8.4 


arc segments to: 

Meshed with elements that span at most maxsegdeg degrees per element 
Boundary property "propname" 

hide: 0 = not hidden in post-processor, 1 == hidden in post processor 

A member of group number group 


A member of the conductor specified by the string "inconductor". If the segment is not part 
of a conductor, this parameter can be specified as "<None>". 


hi_setgroup(n) Set the group associated of the selected items to n. 


Problem Commands 


hi_probdef(units,type,precision,(depth),(minangle)) changes the problem definition. 
The units parameter specifies the units used for measuring length in the problem 
domain. Valid "units" entries are "inches", "millimeters", "centimeters", "mils", "meters, 
and "micrometers". Set problem type to "planar" for a 2-D planar problem, or to "axi" 
for an axisymmetric problem. The precision parameter dictates the precision required 
by the solver. For example, entering 1.E-8 requires the RMS of the residual to be less 
than 10-8. A fourth parameter, representing the depth of the problem in the into-the- 
page direction for 2-D planar problems, can also be specified for planar problems. A sixth 


parameter represents the minimum angle constraint sent to the mesh generator. 


hi_analyze(flag) runs hsolv to solve the problem. The flag parameter controls whether the 
solver window is visible or minimized. For a visible window, specify 0. For a minimized 
window, flag should be set to 1. If no value is specified for flag, the visibility of the solver is 
inherited from the main window, i.e. if the main window is minimized, the solver runs 
minimized, too. 


hi_loadsolution() loads and displays the solution corresponding to the current geometry. 


hi_setfocus(“documentname’’) Switches the heat flow input file upon which Lua 
commands are to act. If more than one heat flow input file is being edited at a time, this 
command can be used to switch between files so that the multiple files can be operated 
upon programmatically via Lua. documentname should contain the name of the desired 
document as it appears on the window’s title bar. 


hi_saveas("filename") saves the file with name "filename". Note if you use a path you must 
use two backslashes e.g. hi_saveas(‘c:\\temp\\myfile.feh’) 
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8.5 


8.6 


Mesh Commands 


hi_createmesh() runs triangle to create a mesh. Note that this is not a necessary precursor of 
performing an analysis, as hi_analyze() will make sure the mesh is up to date before running an 
analysis. Returns number of elements in the mesh. 


hi_showmesh() toggles the flag that shows or hides the mesh. 


hi_purgemesh() clears the mesh out of both the screen and memory. 


Editing Commands 
hi_copyrotate(bx, by, angle, copies, (editaction) ) bx, by — base point for rotation angle — 
angle by which the selected objects are incrementally shifted to make each copy. angle 
is measured in degrees. copies — number of copies to be produced from the selected 


objects. 


editaction 0 —nodes, 1- lines (segments), 2 -block labels, 3 — arc segments, 4- group 


hi_copytranslate(dx, dy, copies, (editaction)) dx,dy — distance by which the selected 
objects are incrementally shifted. copies — number of copies to be produced from 
the selected objects. editaction 0 —nodes, 1 — lines (segments), 2 —block labels, 3 — 


arc segments, 4- group 


e hi_createradius(x,y,r) turns a corner located at (x,y) into a curve of radius r. 


e hi_moverotate(bx,by,shiftangle (,editaction)) 


o bx, by — base point for rotation 
o shiftangle — angle in degrees by which the selected objects are rotated. 
o editaction 

= O-nodes, 

=" 1-lines (segments), 

= 2 —block labels, 

= 3 — arc segments, 


= 4- group 
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8.7 


8.8 


hi_movetranslate(dx,dy,(editaction)) dx,dy — distance by which the selected objects are 
shifted. editaction 0 —nodes, 1 — lines (segments), 2 —block labels, 3 — arc segments, 4- 
group 

hi_scale(bx,by,scalefactor,(editaction)) bx, by — base point for scaling 


scalefactor — a multiplier that determines how much the selected objects are scaled editaction 


0 —-nodes, 1 — lines (segments), 2 —block labels, 3 — arc segments, 4- group 


hi_mirror(x1,y1,x2,y2,editaction)) mirror the selected objects about a line passing 
through the points (x1,y1) and (x2,y2). Valid editaction entries are 0 for nodes, 1 for lines 
(segments), 2 for block labels, 3 for arc segments, and 4 for groups. 


hi_seteditmode(editmode) Sets the current editmode to: 
"nodes" — nodes 

"segments" - line segments 

"arcsegments" - arc segments 

"blocks" - block labels 


"group" - selected group 


This command will affect all subsequent uses of the other editing commands, if they are used 
WITHOUT the editaction parameter. 


Zoom Commands 


|” 


hi_zoomnatural() zooms to a “natural” view with sensible extents. 
hi_zoomout() zooms out by a factor of 50%. 
hi_zoomin() zoom in by a factor of 200%. 


hi_zoom(x1,y1,x2,y2) Set the display area to be from the bottom left corner specified by (x1,y1) to 
the top right corner specified by (x2,y2). 


View Commands 
hi_showgrid() Show the grid points. 
hi_hidegrid() Hide the grid points points. 
hi_gridsnap("flag") Setting flag to “on” turns on snap to grid, setting flag to “off” turns off snap to 


grid. 
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hi_setgrid(density,"type") Change the grid spacing. The density parameter specifies the space 
between grid points, and the type parameter is set to "cart" for Cartesian coordinates or "polar" 
for polar coordinates. 


e hi_refreshview() Redraws the current view. 

e hi_minimize() minimizes the active heat flow input view. 

e hi_maximize() maximizes the active heat flow input view. 

e hi_restore() restores the active heat flow input view from a minimized or maximized state. 


e hi_resize(width,height) resizes the active heat flow input window client area to width x height. 


8.9 Object Properties 


e hi_getmaterial(’materialname’) fetches the material specified by materialname from the 
materials library. 


e hi_addmaterial("materialname", kx, ky, qv, kt) adds a new material with called 
"materialname" with the material properties: kx Thermal conductivity in the x- or r- 
direction. ky Thermal conductivity in the y- or z-direction. qv Volume heat generation 
density in units of W/m? kt Volumetric heat capacity in units of MJ/(m?*K) 


e hi_addpointprop("pointpropname",Tp,qp) adds a new point property of name 
"pointpropname" with either a specified temperature Tp or a point heat generation density 
qp in units of W/m. 


e hi_addboundprop("boundpropname", BdryFormat, Tset, qs, Tinf, h, beta) adds a new 
boundary property with name "boundpropname". 


— For a “Fixed Temperature” type boundary condition, set the Tset parameter to the 
desired temperature and all other parameters to zero. 


— To obtain a “Heat Flux” type boundary condition, set qs to be the heat flux density 
and BdryFormat to 1. Set all other parameters to zero. 


— To obtain a convection boundary condition, set h to the desired heat transfer 
coefficient and Tinf to the desired external temperature and set BdryFormat to 2. 


— Fora Radiation boundary condition, set beta equal to the desired emissivity and Tinf 
to the desired external temperature and set BdryFormat to 3. 


— Fora “Periodic” boundary condition, set BdryFormat to 4 and set all other parameters 
to zero. 
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— For an “Anti-Perodic” boundary condition, set BdryFormat to 5 set all other 
parameters to zero. 


hi_addconductorprop("conductorname", Tc, qc, conductortype) adds a new conductor 
property with name "conductorname" with either a prescribed temperature or a prescribed 
total heat flux. Set the unused property to zero. The conductortype parameter is O for 
prescribed heat flux and 1 for prescribed temperature. 


hi_deletematerial("materialname") deletes the material named "materialname". 


hi_deleteboundprop("boundpropname") deletes the boundary property named 
"boundpropname". 


hi_deleteconductor("conductorname") deletes the conductor named conductorname. 


hi_deletepointprop("pointpropname") deletes the point property named 
"pointpropname" 


hi_modifymaterial("BlockName",propnum,value) This function allows for modification of a 
material’s properties without redefining the entire material (e.g. so that current can be 
modified from run to run). The material to be modified is specified by "BlockName". The 
next parameter is the number of the property to be set. The last number is the value to be 
applied to the specified property. The various properties that can be modified are listed 
below: 


propnum Symbol Description 





BlockName Name of the material 

kx x (or r) direction thermal conductivity 
ky y (or z) direction thermal conductivity 
qs Volume heat generation 


>U NBEO 


kt Volumetric heat capacity 


hi_modifyboundprop("BdryName",propnum,value) This function allows for modification of 
a boundary property. The BC to be modified is specified by "BdryName". The next 
parameter is the number of the property to be set. The last number is the value to be 
applied to the specified property. The various properties that can be modified are listed 
below: 


propnum Symbol Description 





0 BdryName Name of boundary property 
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1 BdryFormat Type of boundary condition: 


O = Prescribed temperature 
1 = Heat Flux 
2 = Convection 
3 = Radiation 
4 = Periodic 
5 = Antiperiodic 
2 Tset Fixed Temperature 
3 qs Prescribed heat flux density 
4 Tinf External temperature 
5 h Heat transfer coefficient 
6 beta Emissivity 


e hi_modifypointprop("PointName",propnum,value) This function allows for modification of 
a point property. The point property to be modified is specified by "PointName". The next 
parameter is the number of the property to be set. The last number is the value to be 
applied to the specified property. The various properties that can be modified are listed 
below: 





propnum Symbol Description 

0 PointName Name of the point property 
1 Tp Prescribed nodal temperature 

2 qp Point heat generation in W/m 


e hi_modifyconductorprop("ConductorName",propnum,value) This function allows for 
modification of a conductor property. The conductor property to be modified is specified 
by "ConductorName". The next parameter is the number of the property to be set. The last 
number is the value to be applied to the specified property. The various properties that can 
be modified are listed below: 





propnum Symbol Description 
0 ConductorName Name of the conductor property 
1 Tc Conductor temperature 
2 qc Total conductor heat flux 
3 ConductorType 0 = Prescribed heat flow, 1 = Prescribed temperature 


e hi_addtkpoint("materialname",T,k) adds the point (T,k) to the thermal conductivity vs. 
temperature curve for the material specified by "materialname". 


e hi_cleartkpoints("materialname") erases all the thermal conductivity points that have been 
defined for the material named "materialname". 
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8.10 Miscellaneous 


hi_savebitmap("filename") saves a bitmapped screenshot of the current view to the file 
specified by "filename”. 


hi_savemetafile("filename") saves a metafile screenshot of the current view to the file 
specified by "filename”. 


hi_refreshview() Redraws the current view. 
hi_close() closes the preprocessor window and destroys the current document. 


hi_shownames(flag) This function allows the user to display or hide the block label names on 
screen. To hide the block label names, flag should be 0. To display the names, the parameter 
should be set to 1. 


hi_readdxf("filename") This function imports a dxf file specified by "filename". 


hi_savedxf("filename") This function saves geometry information in a dxf file specified by 
"filename". 


hi_defineouterspace(Zo,Ro,Ri) defines an axisymmetric external region to be used in 
conjuction with the Kelvin Transformation method of modeling unbounded problems. The 
Zo parameter is the z-location of the origin of the outer region, the Ro parameter is the radius 
of the outer region, and the Ri parameter is the radius of the inner region (i.e. the region of 
interest). In the exterior region, the permeability varies as a function of distance from the 
origin of the external region. These parameters are necessary to define the permeability 
variation in the external region. 


hi_attachouterspace() marks all selected block labels as members of the external region used 
for modeling unbounded axisymmetric problems via the Kelvin Transformation. 


hi_detachouterspace() undefines all selected block labels as members of the external region 
used for modeling unbounded axisymmetric problems via the Kelvin Transformation. 


hi_attachdefault() marks the selected block label as the default block label. This block label 
is applied to any region that has not been explicitly labeled. 


hi_detachdefault() undefines the default attribute for the selected block labels. 


hi_makeABC(n,R,x,y,bc) creates a series of circular shells that emulate the impedance of an 
unbounded domain (i.e. an Improvised Asymptotic Boundary Condition). The n parameter 
contains the number of shells to be used (should be between 1 and 10), R is the radius of the 
solution domain, and (x,y) denotes the center of the solution domain. The bc parameter 
should be specified as O for a Dirichlet outer edge or 1 for a Neumann outer edge. If the 
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function is called without all the parameters, the function makes up reasonable values for 
the missing parameters. 


9 Heat Flow Post Processor Command Set 


There are several Lua scripting commands designed to operate in the postprocessor. As with 
the preprocessor commands, these commands can be used with either the underscore naming 
or with the no-underscore naming convention. 


9.1 Data Extraction Commands 


e ho_lineintegral(type) Calculate the line integral for the defined contour type Integral 





NFO 


3 


Temperature difference (G-t) 

Heat flux through the contour (F :n) 
Contour length 

Average Temperature 


This integral returns either 1 or 2 values, depending on the integral type. 


e ho_blockintegral(type) Calculate a block integral for the selected blocks type Integral 





NFO 
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Average T over the block 

Block Cross-section 

Block Volume 

Average F over the block 4 Average G over the block 


Returns one or two floating point values as results, e.g.: 


Gx, Gy = ho_blockintegral(4) 


e ho_getpointvalues(x,y) Get the values associated with the point at x,y The return values, 
in order, are: 


Symbol Definition 





V 
Fx 


Fy 
Gx 


Gy 
kx 


Temperature 

x- or r- direction component of heat flux density 

y- or z- direction component of heat flux density 

x- or r- direction component of temperature gradient 


y- or z- direction component of temperature gradient 


x- or r- direction component of thermal conductivity 
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ky y- or z- direction component of thermal 
conductivity 


ho_makeplot(PlotType,NumPoints,Filename,FileFormat) Allows Lua access to the X-Y 
plot routines. If only PlotType or only PlotType and NumPoints are specified, the 
command is interpreted as a request to plot the requested plot type to the screen. If, in 
addition, the Filename parameter is specified, the plot is instead written to disk to the 
specified file name as an extended metafile. If the FileFormat parameter is also, the 
command is instead interpreted as a command to write the data to disk to the specfied 
file name, rather than display it to make a graphical plot. Valid entries for PlotType are: 


PlotType Definition 





V (Temperature) 

|D| (Magnitude of heat flux density) 
D.n (Normal heat flux density) 
D.t (Tangential heat flux density) 
|E| (Magnitude of field intensity) 


E . n (Normal field intensity) 


nu BP WN FP OO 


E . t (Tangential field intensity) 


Valid file formats are: 


FileFormat Definition 





Multi-column text with legend 
Multi-column text with no legend 
Mathematica-style formatting 


For example, if one wanted to plot V to the screen with 200 points evaluated to make 
the graph, the command would be: ho_makeplot(0,200) 

If this plot were to be written to disk as a metafile, the command would be: 
ho_makeplot(0,200,"c:\\temp\\temp.emf") 

To write data instead of a plot to disk, the command would be of the form: 
ho_makeplot(0,200,"c:\\temp\\temp.txt",0) 


ho_getprobleminfo() Returns info on problem description. Returns three values: the 
Problem type (0 for planar and 1 for axisymmetric); the depth assumed for planar 
problems in units of meters; and the length unit used to draw the problem, represented 
in meters. 
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9.2 


9.3 


ho_getconductorproperties("conductor") Properties are returned for the conductor 
property named “conductor”. Two values are returned: The temperature of the specified 
conductor, and the total heat flux through the specified conductor. 


Selection Commands 


ho_seteditmode(mode) Sets the mode of the postprocessor to point, contour, or area mode. 


Woot 


Valid entries for mode are "point", "contour", and "area". 
ho_selectblock(x,y) Select the block that contains point (x,y). 


ho_groupselectblock(n) Selects all of the blocks that are labeled by block labels that are members 
of group n. If no number is specified (i.e. ho_groupselectblock()), all blocks are selected. 


ho_selectconductor("name") Selects all nodes, segments, and arc segments that are part of the 
conductor specified by the string ("name"). This command is used to select conductors for the 
purposes of the “weighted stress tensor” force and torque integrals, where the conductors are 
points or surfaces, rather than regions (i.e. can’t be selected with ho_selectblock). 


ho_addcontour(x,y) Adds a contour point at (x,y). If this is the first point then it starts a contour, 
if there are existing points the contour runs from the previous point to this point. The 
ho_addcontour command has the same functionality as a right-button-click contour point 
addition when the program is running in interactive mode. 


ho_bendcontour(angle,anglestep) Replaces the straight line formed by the last two points in 
the contour by an arc that spans angle degrees. The arc is composed of many straight lines, each 
of which is constrained to span no more than anglestep degrees. The angle parameter can take 
on values from -180 to 180 degrees. The anglestep parameter must be greater than zero. If 
there are less than two points defined in the contour, this command is ignored. 


ho_selectpoint(x,y) Adds a contour point at the closest input point to (x,y). If the selected point 
and a previously selected point lie at the ends of an arcsegment, a contour is added that traces 
along the arcsegment. The selectpoint command has the same functionality as the left-button- 
click contour point selection when the program is running in interactive mode. 


ho_clearcontour() Clear a previously defined contour 


ho_clearblock() Clear block selection 


Zoom Commands 


ho_zoomnatural() Zoom to the natural boundaries of the geometry. 


ho_zoomin() Zoom in one level. 
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9.4 


ho_zoomout() Zoom out one level. 


ho_zoom(x1,y1,x2,y2) Zoom to the window defined by lower left corner (x1,y1) and upper right 
corner (x2,y2). 


View Commands 
ho_showmesh() Show the mesh. 
ho_hidemesh() Hide the mesh. 
ho_showpoints() Show the node points from the input geometry. 
ho_hidepoints() Hide the node points from the input geometry. 


ho_smooth("flag") This function controls whether smoothing is applied to the F and G fields which 
are naturally piece-wise constant over each element. Setting flag equal to "on" turns on smoothing 
and setting flag to "off" turns off smoothing. 


ho_showgrid() Show the grid points. 


ho_hidegrid() Hide the grid points. 
ho_gridsnap("flag") Setting flag to “on” turns on snap to grid, setting flag to “off” turns off snap 


to grid. 


ho_setgrid(density,"type") Change the grid spacing. The density parameter specifies the space 
between grid points, and the type parameter is set to "cart" for Cartesian coordinates or "polar" 
for polar coordinates. 


ho_hidedensityplot() hides the heat flux density plot. 


ho_showdensityplot(legend,gscale,type,upper,lower) Shows the heat flux density plot with 
options: legend Set to 0 to hide the plot legend or 1 to show the plot legend. gscale Set to O for a 
colour density plot or 1 for a grey scale density plot. upper Sets the upper display limit for the 
density plot. lower Sets the lower display limit for the density plot. 

type Sets the type of density plot. A value of O plots temperature, 1 plots the magnitude of F, 
and 2 plots the magnitude of G 


ho_hidecontourplot() Hides the contour plot. 


ho_showcontourplot(numcontours,lower V,upper V) shows the V contour plot with options: 
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9.5 


numcontours Number of equipotential lines to be plotted. 
upper V Upper limit for contours. lower V Lower limit for 
contours. 


If ho_numcontours is -1 all parameters are ignored and default values are used, e.g. 
show contour plot(-1) 


ho_showvectorplot(type,scalefactor) controls the display of vectors denoting the field strength 
and direction. The parameters taken are the type of plot, which should be set to O for no vector 
plot, 1 for heat flux density F, and 2 for temperature gradient G. The scalefactor determines the 
relative length of the vectors. If the scale is set to 1, the length of the vectors is chosen so that the 
highest flux density corresponds to a vector that is the same length as the current grid size setting. 


ho_minimize() minimizes the active heat flow input view. 
ho_maximize() maximizes the active heat flow input view. 
ho_restore() restores the active heat flow input view from a minimized or maximized state. 


ho_resize(width,height) resizes the active heat flow input window client area to width x height. 


Miscellaneous 


ho_close() close the current postprocessor window. 
ho_refreshview() Redraws the current view. 
ho_reload() Reloads the solution from disk. 


ho_savebitmap("filename") saves a bitmapped screen shot of the current view to the file specified 
by "filename". 


ho_savemetafile("filename") saves a metafile screenshot of the current view to the file specified 
by "filename”. 


ho_shownames(flag) This function allows the user to display or hide the block label names on 
screen. To hide the block label names, flag should be 0. To display the names, the parameter should 
be set to 1. 


ho_numnodes Returns the number of nodes in the in-focus heat flow output mesh. 
ho_numelements Returns the number of elements in the in-focus heat flow output mesh. 
ho_getnode(n) Returns the (x,y) or (r,z) position of the nth mesh node. 

ho_getelement(n) returns the following properties for the nth element: 
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10 


. Index of first element node 

. Index of second element node 

. Index of third element node 

. x (or r) coordinate of the element centroid 
. y (or z) coordinate of the element centroid 


. element area using the length unit defined for the problem 


N DO UW BP WN PR 


. group number associated with the element 


Current Flow Preprocessor Lua Command Set 


Several different commands are available in the preprocessor. Two naming conventions can 
be used: one which separates words in the command names by underscores, and one that 
eliminates the underscores. 


10.1 


10.2 


Object Add/Remove Commands 


ci_addnode(x,y) Add a new node at x,y 


ci_addsegment(x1,y1,x2,y2) Add a new line segment from node closest to (x1,y1) to node closest 
to (x2,y2) 


ci_addblocklabel(x,y) Add a new block label at (x,y) 


ci_addarc(x1,y1,x2,y2,angle, maxseg) Add a new arc segment from the nearest node to (x1,y1) to 
the nearest node to (x2,y2) with angle ‘angle’ divided into ‘maxseg’ segments. 


ci_deleteselected Delete all selected objects. 
ci_deleteselectednodes Delete selected nodes. 
ci_deleteselectedlabels Delete selected block labels. 
ci_deleteselectedsegments Delete selected segments. 


ci_deleteselectedarcsegments Delete selects arcs. 


Geometry Selection Commands 


ci_clearselected() Clear all selected nodes, blocks, segments and arc segments. 
ci_selectsegment(x,y) Select the line segment closest to (x,y) 


ci_selectnode(x,y) Select the node closest to (x,y). Returns the coordinates of the selected node. 
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e ci_selectlabel(x,y) Select the label closet to (x,y). Returns the coordinates of the selected label. 


e ci_selectarcsegment(x,y) Select the arc segment closest to (x,y) 


e ci_selectgroup(n) Select the n“ group of nodes, segments, arc segments and block labels. 
This function will clear all previously selected elements and leave the edit mode in 4 (group) 


e ci_selectcircle(x,y,R,editmode) selects objects within a circle of radius R centered at (x,y). If only 
x, y, and R paramters are given, the current edit mode is used. If the editmode parameter is 
used, 0 denotes nodes, 2 denotes block labels, 2 denotes segments, 3 denotes arcs, and 4 
specifies that all entity types are to be selected. 


e ci_selectrectangle(x1,y1,x2,y2,editmode) selects objects within a rectangle defined by points 
(x1,y1) and (x2,y2). If no editmode parameter is supplied, the current edit mode is used. If the 
editmode parameter is used, 0 denotes nodes, 2 denotes block labels, 2 denotes segments, 3 
denotes arcs, and 4 specifies that all entity types are to be selected. 


10.3 Object Labeling Commands 


e ci_setnodeprop("propname",groupno, "inconductor") Set the selected nodes to have the 
nodal property "propname" and group number groupno. The "inconductor" string specifies 
which conductor the node belongs to. If the node doesn’t belong to a named conductor, this 
parameter can be set to "<None>". 


e ci_setblockprop("blockname", automesh, meshsize, group) Set the selected block labels to 


have the properties: Block property "blockname". 


automesh: 0 = mesher defers to mesh size constraint defined in meshsize, 1 = mesher 
automatically chooses the mesh density. meshsize: size constraint on the mesh in the 


block marked by this label. A member of group number group 
e ci_setsegmentprop("propname", elementsize, automesh, hide, group, "inconductor",) Set 
the select segments to have: 
Boundary property "propname" 
Local element size along segment no greater than elementsize 


automesh: 0 = mesher defers to the element constraint defined by elementsize, 1 = 
mesher automatically chooses mesh size along the selected segments hide: 0 = not 
hidden in post-processor, 1 == hidden in post processor 


A member of group number group 


A member of the conductor specified by the string "inconductor". If the segment is not part 
of a conductor, this parameter can be specified as "<None>". 
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e ci_setarcsegmentprop(maxsegdeg, "propname", hide, group, "inconductor") Set the 
selected arc segments to: 


Meshed with elements that span at most maxsegdeg degrees per element 
Boundary property "propname" 

hide: 0 = not hidden in post-processor, 1 == hidden in post processor 

A member of group number group 


A member of the conductor specified by the string "inconductor". If the segment is not part 
of a conductor, this parameter can be specified as "<None>". 


e ci_setgroup(n) Set the group associated of the selected items to n. 


10.4 Problem Commands 


e ci_probdef(units,type,frequency,precision,(depth),(minangle)) changes the problem 
definition. The units parameter specifies the units used for measuring length in the problem 
domain. Valid "units" entries are "inches", "millimeters", "centimeters", "mils", "meters, 
and "micrometers". Set problemtype to "planar" for a 2-D planar problem, or to "axi" for 
an axisymmetric problem. The frequency parameter specifies the frequency in Hz at which 
the analysis ois to be performed. The precision parameter dictates the precision required 
by the solver. For example, entering 1.E-8 requires the RMS of the residual to be less than 
10-8. A fourth parameter, representing the depth of the problem in the into-the-page 
direction for 2-D planar problems, can also be specified for planar problems. A sixth 


parameter represents the minimum angle constraint sent to the mesh generator. 


e ci_analyze(flag) runs belasolv to solve the problem. The flag parameter controls whether 
the solver window is visible or minimized. For a visible window, specify 0. For a minimized 
window, flag should be set to 1. If no value is specified for flag, the visibility of the solver is 
inherited from the main window, i.e. if the main window is minimized, the solver runs 
minimized, too. 


e ci_loadsolution() loads and displays the solution corresponding to the current geometry. 


e ci_setfocus("documentname") Switches the electrostatics input file upon which Lua 
commands are to act. If more than one electrostatics input file is being edited at a time, 
this command can be used to switch between files so that the mutiple files can be operated 
upon programmatically via Lua. documentname should contain the name of the desired 
document as it appears on the window’s title bar. 


e ci_saveas("filename") saves the file with name "filename". Note if you use a path you must 
use two backslashes e.g. ci_saveas(‘c:\\temp\\myfemmfile.fee’) 
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10.5 


10.6 


Mesh Commands 


ci_createmesh() runs triangle to create a mesh. Note that this is not a necessary precursor of 
performing an analysis, as ci_analyze() will make sure the mesh is up to date before running an 
analysis. The number of elements in the mesh is pushed back onto the lua stack. 


ci_showmesh() toggles the flag that shows or hides the mesh. 


ci_purgemesh() clears the mesh out of both the screen and memory. 


Editing Commands 


ci_copyrotate(bx, by, angle, copies, (editaction) ) bx, by — base point for rotation angle 
— angle by which the selected objects are incrementally shifted to make each copy. 


angle is measured in degrees. 


copies — number of copies to be produced from the selected objects. editaction 0 — 


nodes, 1 —lines (segments), 2 -block labels, 3 — arc segments, 4- group 


ci_copytranslate(dx, dy, copies, (editaction)) 
dx,dy — distance by which the selected objects are incrementally shifted. copies — 
number of copies to be produced from the selected objects. editaction 0 —nodes, 


1- lines (segments), 2 -block labels, 3 — arc segments, 4- group 
ci_createradius(x,y,r) turns a corner located at (x,y) into a curve of radius r. 


ci_moverotate(bx,by,shiftangle (editaction)) bx, by — base point for rotation shiftangle 
— angle in degrees by which the selected objects are rotated. editaction 0 —nodes, 1- 


lines (segments), 2 -block labels, 3 — arc segments, 4- group 
ci_movetranslate(dx,dy,(editaction)) dx,dy — distance by which the selected objects are 
shifted. editaction 0 —nodes, 1 — lines (segments), 2 —block labels, 3 — arc segments, 4- 
group 

ci_scale(bx,by,scalefactor,(editaction)) bx, by — base point for scaling 


scalefactor — a multiplier that determines how much the selected objects are scaled editaction 


0 —-nodes, 1 — lines (segments), 2 —block labels, 3 — arc segments, 4- group 


54 


10.7 


10.8 


ci_mirror(x1,y1,x2,y2,(editaction)) mirror the selected objects about a line passing 
through the points (x1,y1) and (x2,y2). Valid editaction entries are O for nodes, 1 for 
lines (segments), 2 for block labels, 3 for arc segments, and 4 for groups. 


ci_seteditmode(editmode) Sets the current editmode to: 
"nodes" — nodes 

"segments" - line segments 

“arcsegments" - arc segments 

"blocks" - block labels 


"group" - selected group 


This command will affect all subsequent uses of the other editing commands, if they are used 
WITHOUT the editaction parameter. 


Zoom Commands 
ci_zoomnatural() zooms to a “natural” view with sensible extents. 
ci_zoomout() zooms out by a factor of 50%. 
ci_zoomin() zoom in by a factor of 200%. 
ci_zoom(x1,y1,x2,y2) Set the display area to be from the bottom left corner specified by (x1,y1) to 


the top right corner specified by (x2,y2). 


View Commands 


ci_showgrid() Show the grid points. 
ci_hidegrid() Hide the grid points points. 


ci_gridsnap("flag") Setting flag to ”on” turns on snap to grid, setting flag to ”off’turns off snap to 
grid. 


ci_setgrid(density,"type") Change the grid spacing. The density parameter specifies the space 
between grid points, and the type parameter is set to "cart" for Cartesian coordinates or "polar" 
for polar coordinates. 


ci_refreshview() Redraws the current view. 
ci_minimize() minimizes the active magnetics input view. 


ci_maximize() maximizes the active magnetics input view. 
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e ci_restore() restores the active magnetics input view from a minimized or maximized state. 


e ci_resize(width,height) resizes the active magnetics input window client area to width x height. 


10.9 Object Properties 


e ci_getmaterial(’materialname’) fetches the material specified by materialname from the 
materials library. 
e ci_addmaterial("materialname", ox, oy, ex, ey, Itx, Ity) adds a new material with called 


"materialname" with the material properties: 


ox Electrical conductivity in the x- or r-direction in units of S/m. 
oy Electrical conductivity in the y- or z-direction in units of S/m. 
ex Relative permittivity in the x- or r-direction. 
ey Relative permittivity in the y- or z-direction. 
Itx Dielectric loss tangent in the x- or r-direction. 
Ity Dielectric loss tangent in the y- or z-direction. 
e ci_addpointprop("pointpropname",Vp,jp) adds a new point property of name 


"pointpropname" with either a specified potential Vp a point current density jp in units of 
A/m. 


e ci_addboundprop("boundpropname", Vs, js, cO, c1, BdryFormat) adds a new boundary 
property with name "boundpropname" 


For a “Fixed Voltage” type boundary condition, set the Vs parameter to the desired voltage 
and all other parameters to zero. 


To obtain a “Mixed” type boundary condition, set C1 and CO as required and BdryFormat to 
1. Set all other parameters to zero. 


To obtain a prescribes surface current density, set js to the desired current density in A/m? 
and set BdryFormat to 2. 


For a “Periodic” boundary condition, set BdryFormat to 3 and set all other parameters to 
zero. 


For an “Anti-Perodic” boundary condition, set BdryFormat to 4 set all other parameters to 
zero. 


e ci_addconductorprop("conductorname", Vc, jc, conductortype) adds a new conductor 
property with name "conductorname" with either a prescribed voltage or a prescribed 
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total current. Set the unused property to zero. The conductortype parameter is O for 
prescribed current and 1 for prescribed voltage. 


ci_deletematerial("materialname") deletes the material named "materialname". 


ci_deleteboundprop("boundpropname") deletes the boundary property named 
"boundpropname". 


ci_deleteconductor("conductorname") deletes the conductor named conductorname. 


ci_deletepointprop("pointpropname") deletes the point property named 
"pointpropname" 


ci_modifymaterial("BlockName",propnum,value) This function allows for modification 
of a material’s properties without redefining the entire material (e.g. so that current 
can be modified from run to run). The material to be modified is specified by 
"BlockName". The next parameter is the number of the property to be set. The last 
number is the value to be applied to the specified property. The various properties that 
can be modified are listed below: 


propnum Symbol Description 





BlockName Name of the material 

ox x (or r) direction conductivity 

oy y (or z) direction conductivity 

ex x (or r) direction relative permittivity 
ey y (or z) direction relative permittivity 
Itx x (or r) direction dielectric loss tangent 
Ity y (or z) direction dielectric loss tangent 
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ci_modifyboundprop("BdryName",propnum,value) This function allows for 
modification of a boundary property. The BC to be modified is specified by "BdryName". 
The next parameter is the number of the property to be set. The last number is the 
value to be applied to the specified property. The various properties that can be 
modified are listed below: 





propnum Symbol Description 

0 BdryName Name of boundary property 
1 Vs Fixed Voltage 

2 js Prescribed current density 
3 c0 Mixed BC parameter 

4 c1 Mixed BC parameter 
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5 BdryFormat Type of boundary condition: 
O = Prescribed V 


1 = Mixed 

2 = Surface current density 
3 = Periodic 

4 = Antiperiodic 


e ci_modifypointprop("PointName",propnum,value) This function allows for 
modification of a point property. The point property to be modified is specified by 
"PointName". The next parameter is the number of the property to be set. The last 
number is the value to be applied to the specified property. The various properties that 
can be modified are listed below: 





propnum Symbol Description 
0 PointName Name of the point property 
1 Vp Prescribed nodal voltage 
2 jp Point current density in A/m 


e ci_modifyconductorprop("ConductorName",propnum,value) This function allows for 
modification of a conductor property. The conductor property to be modified is 
specified by "ConductorName". The next parameter is the number of the property to 
be set. The last number is the value to be applied to the specified property. The various 
properties that can be modified are listed below: 





propnum Symbol Description 

0 ConductorName Name of the conductor property 

1 Vc Conductor voltage 

2 jc Total conductor current 

3 ConductorType 0 = Prescribed current, 1 = Prescribed voltage 


10.10 Miscellaneous 


e ci_savebitmap("filename") saves a bitmapped screenshot of the current view to the file 
specified by "filename”. 


e ci_savemetafile("filename") saves a metafile screenshot of the current view to the file 
specified by "filename”. 


e ci_refreshview() Redraws the current view. 


e ci_close() closes the preprocessor window and destroys the current document. 
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ci_shownames(flag) This function allows the user to display or hide the block label names on 
screen. To hide the block label names, flag should be 0. To display the names, the parameter 
should be set to 1. 


ci_readdxf("filename") This function imports a dxf file specified by "filename". 


ci_savedxf("filename") This function saves geometry information in a dxf file specified by 
"filename". 


ci_defineouterspace(Zo,Ro,Ri) defines an axisymmetric external region to be used in 
conjunction with the Kelvin Transformation method of modeling unbounded problems. The 
Zo parameter is the z-location of the origin of the outer region, the Ro parameter is the 
radius of the outer region, and the Ri parameter is the radius of the inner region (i.e. the 
region of interest). In the exterior region, the permeability varies as a function of distance 
from the origin of the external region. These parameters are necessary to define the 
permeability variation in the external region. 


ci_attachouterspace() marks all selected block labels as members of the external region used 
for modeling unbounded axisymmetric problems via the Kelvin Transformation. 


ci_detachouterspace() undefines all selected block labels as members of the external region 
used for modeling unbounded axisymmetric problems via the Kelvin Transformation. 


ci_attachdefault() marks the selected block label as the default block label. This block label 
is applied to any region that has not been explicitly labeled. 


ci_detachdefault() undefines the default attribute for the selected block labels. 


ci_makeABC(n,R,x,y,bc) creates a series of circular shells that emulate the impedance of an 
unbounded domain (i.e. an Improvised Asymptotic Boundary Condition). The n parameter 
contains the number of shells to be used (should be between 1 and 10), R is the radius of 
the solution domain, and (x,y) denotes the center of the solution domain. The bc parameter 
should be specified as O for a Dirichlet outer edge or 1 for a Neumann outer edge. If the 
function is called without all the parameters, the function makes up reasonable values for 
the missing parameters. 


Current Flow Post Processor Command Set 


There are several Lua scripting commands designed to operate in the postprocessor. As with 
the preprocessor commands, these commands can be used with either the underscore 
naming or with the no-underscore naming convention. 
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11.1 Data Extraction Commands 


e co_lineintegral(type) Calculate the line integral for the defined contour type 


Integral 
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E-t 

J:n 

Contour length 

Average voltage over contour 
Force from stress tensor 
Torque from stress tensor 


This integral return either 1 or 2 values, depending on the integral type 


e co_blockintegral(type) Calculate a block integral for the selected blocks type 


Integral 
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m. 
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Real Power 

Reactive Power 

Apparent Power 

Time-Average Stored Energy 

Block cross-section area 

Block volume 

x (or r) direction Weighted Stress Tensor force, DC component 

y (or z) direction Weighted Stress Tensor force, DC component 

x (or r) direction Weighted Stress Tensor force, 2x frequency component 
y (or z) direction Weighted Stress Tensor force, 2x frequency component 
Weighted Stress Tensor torque, DC component 

Weighted Stress Tensor torque, 2x frequency component Returns a value that can be 


complex, as necessary. 


e co_getpointvalues(x,y) Get the values associated with the point at x,y The 
return values, in order, are: 


Symbol Definition 





V 

Jx 
Jy 
Kx 
Ky 


Voltage 

x- or r- direction component of current density 
y- or z- direction component of current density 
x- or r- direction component of AC conductivity 
y- or z- direction component of AC conductivity 
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Ex x- or r- direction component of electric field intensity 


Ey y- or z- direction component of electric field intensity 

ex x- or r- direction component of permittivity 

ey y- or z- direction component of permittivity 

Jdx x- or r- direction component of displacement current density 
Jdy y- or z- direction component of displacement current density 
Ox x- or r- direction component of permittivity 

oy y- or z- direction component of permittivity 

Jcx x- or r- direction component of conduction current density 
Jcy y- or z- direction component of conduction current density 


¢co_makeplot(PlotType,NumPoints, Filename, FileFormat) Allows Lua access to the X-Y 
plot routines. If only PlotType or only PlotType and NumPoints are specified, the 
command is interpreted as a request to plot the requested plot type to the screen. If, 
in addition, the Filename parameter is specified, the plot is instead written to disk to 
the specified file name as an extended metafile. If the FileFormat parameter is also, the 
command is instead interpreted as a command to write the data to disk to the specfied 
file name, rather than display it to make a graphical plot. Valid entries for PlotType are: 


PlotType Definition 





0 V (Voltage) 

1 |J| (Magnitude of current density) 

2 J . n (Normal current density) 

3 J . t (Tangential current density) 

4 JE] (Magnitude of field intensity) 

5 E . n (Normal field intensity) 

6 E . t (Tangential field intensity) 

7 |Jc| (Magnitude of conduction current density) 
8 Jc . n (Normal conduction current density) 

9 Jc . t (Tangential conduction current density) 
10 |Jd| (Magnitude of displacement current density) 
11 Jd . n (Normal displacement current density) 
12 Jd . t (Tangential displacement current density) 


Valid file formats are: 


FileFormat Definition 





0 Multi-column text with legend 
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Multi-column text with no legend 
Mathematica-style formatting 


For example, if one wanted to plot V to the screen with 200 points evaluated to make the graph, 
the command would be: co_makeplot(0,200) 

If this plot were to be written to disk as a metafile, the command would be: 
co_makeplot(0,200,"c:\\temp\\temp.emf") 

To write data instead of a plot to disk, the command would be of the form: 
co_makeplot(0,200,"c:\\temp.txt\\",0) 


e co_getprobleminfo() Returns info on problem description. Returns four values: 








Return value Definition 
1 problem type 

2 frequency in Hz 

3 depth assumed for planar problems in meters 

4 length unit used to draw the problem, represented in meters. 


e co_getconductorproperties("conductor") Properties are returned for the conductor 
property named “conductor”. Two values are returned: The voltage of the specified 
conductor, and the current on the specified conductor. 


11.2 Selection Commands 


e co_seteditmode(mode) Sets the mode of the postprocessor to point, contour, or area mode. 


wou 


Valid entries for mode are "point", "contour", and "area". 
e co_selectblock(x,y) Select the block that contains point (x,y). 


e co_groupselectblock(n) Selects all of the blocks that are labeled by block labels that are 
members of group n. If no number is specified (i.e. co_groupselectblock()), all blocks are 
selected. 


e co_selectconductor("name") Selects all nodes, segments, and arc segments that are part of 
the conductor specified by the string ("name"). This command is used to select conductors 
for the purposes of the “weighted stress tensor” force and torque integrals, where the 
conductors are points or surfaces, rather than regions (i.e. can’t be selected with 
co_selectblock). 


e co_addcontour(x,y) Adds a contour point at (x,y). If this is the first point then it starts a 
contour, if there are existing points the contour runs from the previous point to this point. 
The co_addcontour command has the same functionality as a right-button-click contour point 
addition when the program is running in interactive mode. 
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e co_bendcontour(angle,anglestep) Replaces the straight line formed by the last two points in 
the contour by an arc that spans angle degrees. The arc is composed of many straight lines, 
each of which is constrained to span no more than anglestep degrees. The angle parameter 
can take on values from -180 to 180 degrees. The anglestep parameter must be greater than 
zero. If there are less than two points defined in the contour, this command is ignored. 


e co_selectpoint(x,y) Adds a contour point at the closest input point to (x,y). If the selected 
point and a previously selected point lie at the ends of an arcsegment, a contour is added 
that traces along the arcsegment. The selectpoint command has the same functionality as 
the left-button-click contour point selection when the program is running in interactive 
mode. 


e co_clearcontour() Clear a prevously defined contour 


e co_clearblock() Clear block selection 


11.3 Zoom Commands 
e co_zoomnatural() Zoom to the natural boundaries of the geometry. 
e co_zoomin() Zoom in one level. 
e coO_zoomout() Zoom out one level. 
e co_zoom(x1,y1,x2,y2) Zoom to the window defined by lower left corner (x1,y1) and upper right 


corner (x2,y2). 


11.4 View Commands 


e co_showmesh() Show the mesh. 

e co_hidemesh() Hide the mesh. 

e co_showpoints() Show the node points from the input geometry. 
e co_hidepoints() Hide the node points from the input geometry. 


e co_smooth("flag") This function controls whether or not smoothing is applied to the D and E fields 
which are naturally piece-wise constant over each element. Setting flag equal to "on" turns on 
smoothing and setting flag to "off" turns off smoothing. 


e co_showgrid() Show the grid points. 


e co_hidegrid() Hide the grid points points. 
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co_gridsnap("flag") Setting flag to “on” turns on snap to grid, setting flag to ”off” turns off 
snap to grid. 


co_setgrid(density,"type") Change the grid spacing. The density parameter specifies the space 
between grid points, and the type parameter is set to "cart" for Cartesian coordinates or "polar" 
for polar coordinates. 


co_hidedensityplot() hides the current density plot. 


co_showdensityplot(legend,gscale,type,upper,lower) Shows the current density plot with options: 
legend Set to 0 to hide the plot legend or 1 to show the plot legend. gscale Set to O for a colour 
density plot or 1 for a grey scale density plot. upper Sets the upper display limit for the density 
plot. lower Sets the lower display limit for the density plot. type Sets the type of density plot. 
Specific choices for the type of density plot include: 


type Description 





[V| 
|Re(V)| 
[Im(V)| 
II 
|Re()| 
Imi) 
[EI 
|Re(E)| 
|Im(E)| 


co_hidecontourplot() Hides the contour plot. 


CON DU BPWN FP OO 


co_showcontourplot(numcontours,lower V,upper V),type shows the V contour plot with 
options: 


numcontours Number of equipotential lines to be 
plotted; upper V Upper limit for contours; lower V 
Lower limit for contours; type the type of contour 
plot to be rendered. 


If co_numcontours is -1 all parameters are ignored and default values are used, 
e.g. show contour plot(-1) 

The type can take on the values of "real", "imag", or "both", denoting the real part of voltage, 
the imaginary part of voltage, or both components of voltage. 
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11.5 


co_showvectorplot(type,scalefactor) controls the display of vectors denoting the field strength 
and direction. The type parameter can take on the following values: 


type Description 





Re(J) and Im(J) 
Re(E) and /m(E) 


0 No vector plot 
1 Re(J) 

2 Re(E) 

3 Im(J) 

4 Im(E) 

5 

6 


The scalefactor determines the relative length of the vectors. If the scale is set to 1, the length of 
the vectors is chosen so that the highest field magnitude corresponds to a vector that is the same 
length as the current grid size setting. 


co_minimize() minimizes the active magnetics input view. 
co_maximize() maximizes the active magnetics input view. 
co_restore() restores the active magnetics input view from a minimized or maximized state. 


co_resize(width,height) resizes the active magnetics input window client area to width x height. 


Miscellaneous 


co_close() close the current postprocessor window. 
co_refreshview() Redraws the current view. 
co_reload() Reloads the solution from disk. 


co_savebitmap("filename") saves a bitmapped screen shot of the current view to the file specified 
by "filename". 


co_savemetafile("filename") saves a metafile screenshot of the current view to the file specified 
by "filename”. 


co_shownames(flag) This function allows the user to display or hide the block label names on 
screen. To hide the block label names, flag should be 0. To display the names, the parameter should 
be set to 1. 


co_numnodes Returns the number of nodes in the in focus current flow output mesh. 
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e co_numelements Returns the number of elements in the in focus current flow output mesh. 


e co_getnode(n) Returns the (x,y) or (r,z) position of the nth mesh node. 


e co_getelement(n) returns the following properties for the nth element: 


1. Index of first element node 

2. Index of second element node 
3. 
4 


. X (or r) coordinate of the element centroid 


Index of third element node 


5. y (or z) coordinate of the element centroid 


. element area using the length unit defined for the problem 


. group number associated with the element 
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